ADA0Z71  40 


DNA  3827T 


COLLISIONAL  MOMENTUM  AND  ENERGY 
TRANSFER  RATES  FOR  TWO-FLUID 
NUCLEAR-BURST  SIMULATIONS 


Mission  Research  Corporation 

735  State  Street 

Santa  Barbara,  California  93101 


25  July  1975 


Topical  Report  for  Period  October  1974— July  1975 


CONTRACT  No.  DNA  001-74-C-0144 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THIS  WORK  SPONSORED  BY  THE  DEFENSE  NUCLEAR 
UNDER  SUBTASK  S99QAXHC06545. 


4 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  D.  C.  20305 


'•  ''-wwswwa. 


UNCLASS  iriLD 

~t  • Jl.',  ’ . '•  >Ai^'  1 1 - A *■»  >N  *■  1 •*!,!>/.'.»  ■H*,fi.  I*n,n  J !,(.'(•  -• 


( / / REPORT  DOCUMENTATION  PAGE 

oj  "i~ ; TTTSguw  ^ CT7  ’ “ '.'  Y.O\7 


DApp  Ki.AI?  INSTRUCTIONS 

.___! UM'OKK  ( uMI'i.l.TINT.  IORM 

. l,OVT  ACCl.'iMOM  l.o  i Mni'inr' J f..ATALl)T  »ltlMlll.n 


m 


COLL  I SIGNAL  MOMENTUM  AND  ENERGY  TRANSTLR  RATES 
FOR  TWO-FLUID  NUCLF.AR-BUkST  SIMULATIONS* 


— ) 

oec 

aga 

b 

K IX 

t ' 

U (I  It  I 'nM'N  % •>»<  .r  A!.'  .*  A T :0»1  N»Mi  AND  A'/(jN  I'  5 ', 

Mission  Research  Corporation 
735  State  Street 

Santa  Sarbara,  California  93101 

’*  f'Pjl  RO»  LINO  O I ’ t N AMI  AND  ADOmsS 

Di  rector 

Defense  Nuclear  Agency 
Washington,  D.C._  20305 

7 amT'a  ad? 


HWOli^AM  [ LlMLNf  PROJECT  TALK 

AI>LA  <»  A'/l-K  UNIT  NjMOEHS 


NY/1SD  Subtask  S95QAXHC065-OS 


25  July:i975 


1-1 

\ 

L 

' IS/i 

- ■ JM"  v <_  l AS.  ' fMs 

UNCLASSIFIED 

DF.'UAWlFK  ATI  ON  DOWNiT.HADlNG 
SCHEDULE 

•G  Dl  jT  i<1  at.  T|ON  S T A T I!  M f . M 

Appr.oved  JFor  pub  I ic  release;  distribution  unlimited. 

//■'-  ) i A’  i )-  A)  ;'iVM  \ 


V6j  i ^2  - Ny!}  ;:k:  Q i' X H _\ 


’7  M'jf  iMilj’  ION  ST  ATEMCN  T '•-/  thr  nhs!tnct  ••fifi'fr  / »i  tilm'k  JO,  tl  ./» / f «»r •*/ if  from  fiv/utrO 

(fn)  C4ATI 


Y s /*  i 


'0  WfBi.E»EHTAItY  H’li 


This  work  sponsored  by  the  Defense  Nuclear  Agency  under  Subtask 
S99QAXUC065-05 . 


'Ey  ADROS  1.  •ri'i/i  . 


i.i  fn,  M ?•/.*.  ^ iiiiwhi'r , 


Atomic  Collisions  Charge-transfer  Cross  Section 

High-altitude  Nuclear  Bursts  Elastic  Scattering  Cross  Section 

Two-fluid  Magnetohydrodynami c Simulations  Inelastic  Scattering  Cross  Section 
Momentum  Coupling  Coefficients  Heat-transfer  Coefficients 

ii(tST  ft  ACT  >n  rtu-iTw  '■tlr  H net  •-•  s .son  twit  At-  M *<  h ; 

""■’A  compendium  of  cross  sections  is  given  for  atomic  collision  processes 

relevant  to  high-altitude  nuclear  bursts.  Experimental  or  theoretical  cross 
sections  for  elastic  scattering,  inelastic  scattering,  and  charge  exchange 
in  ion-neutral  and  neutral -neutral  collisions  and  for  elastic  scattering  in 
electron-neutral  and  electron-ion  collisions  are  given  for  atmospheric 
species  in  the  velocity  range  ofi  10^  to  10^  cm/sec.  The  electrical  resis- 
tivity and  the  rates  of  momentum  and  heat  transfer  between  the  ion  fluid 
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■■and  the  neutral  fluid  are  obtained  from  these  cross  sections  for  the  case 
that  the  particle  velocity  distribution  of  each  species  is  Maxwell- 
Boltzmann.  The  coefficients  for  these  quantities  are  given  graphically 
and  numerically  for  temperatures  in  the  range  0.01  to  500  eV.  A summary 
shows  how  to  use  these  data  in  two-fluid  magnetohydrodynamic  simulations 
of  high-altitude  nuclear  bursts  at  tines  subsequent  to  one  second. 
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IN1RODDCI ION  AND  SUMMARY 


Analysis  of  (In-  effect  ul  high- a 1 1 1 1 tide  nuclear  bursts  on  coiiiuuihi- 
cation,  radar,  optical,  and  infrared  systems  requires  reliable  computations 
of: 

(1)  the  interaction  of  the  hir.h  velocity  weapon  debris  with  the 
surround itn;  air; 

UJ  the  subsequent  motion  of  the  heated  .and  ionir.cd  air; 

the  late- time  t\»riuat  ion.  vl  r i ft  , and  decav  of  striatious  in 
hi  j'.lt-a ! t i t ude  ion-electron  plasma. 

lor  hi  qh-nl  t i t tide  nuclear  bursts  such  phenomena  depend  strong  l > oil  IK  tales 
of  collisions  between  ions  and  neutrals,  electrons  and  neutrals,  ant  ’ -crrotv 
and  ions,  lor  hursts  below  I on  Km  altitude  collisions  are  suff i c i < at  i \ 
rapid  that  the  one- fluid  hydrodynamic  equations  describe  most  pin  noincna 
adequately,  lor  bursts  above  lt»u  Km  altitude,  however,  the  lower  densities 
cause  the  col  lision  mean  tree  paths  to  he  quite  K.ni;.  an.1  so  the  detail'-  of 
the  collision  processes  In-come  important.  ihis  report  presents  a survey  of 
the  collision  cross  sect  ions  relevant  to  nsl.-.i’’.  t <Ji  busts  and  describes 
bow  to  use  this  data  in  two-fluid  iiur.iietohydrodynanuc  i.'lidi  s titulat  ions  for 
prediction  of  burst  phenomena. 


the  relevant  types  of  collisions  arc  elastic  s>-at  tori  up,,  inelastic 
'•cat  t er  i m; , and  chare.e  e.vch.mro.  \ reliable  analysis  of  burst  phenomena 

■I  ■> 

requires  these  cross  sections  mar  the  velocity  ranr.e  of  It'  to  It’  cui,' sec. 

l’uriiu;  the  first  second,  much  of  the  burst  debris  Kinetic  encre.y  is  traits 

'erred  to  the*  snrroimdini;  air,  creat  in;;  air  ions  with  velocit  i es  .is  larre  as 
»> 

- Ilf  ein/sec.  the  subsequent  expansion  and  rise  ot  these  t i reha  I I ions  launch 
a netit  ra  1 - f lu  i d shock  from  the  burst  aria  th.it  involves  particle  velocities 
as  low  as  - It’'  cm/ see.  The  sources  of  those  cross-sect  ion  data  a ail  relevant 
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theoretical  >. « >111 1 >i  1 1 .•  i ><•■■-.  .in-  di".v  r i bed  in  Sections  d,  ' ;hkI  Appendix  C. 
ligurcs  I I through  I S miiiim.tr  i ;".c  the  cross  scU  inn  il.it  a fur  the  panicle 
sivcii's  import  .‘Mir  in  high  alt  truck'  hursts. 

la  use  the  e i’iiss-  sect  mu  dal  a.  one  must  define'  the  phenomena  U>  l>e 
anal.'.ed  and  then  do  the  appropriate  integrals  over  the  particle  velocity 
distributions.  Huring  the  first  second,  is  hen  the  particle  velocities  arc 
extreme!;.  high  ( * lit  cm/seel,  these  phenomena  tend  to  involve  very  special  ir.od 
applis.it  ions  (for  example  (llll.X  losses  and  col  I isiontil  pici.up),  which  ice  do 
not  attempt  to  analy-e  here.  Subsequent  to  one  second  the  main  applieulion 
is  a tivo  fluid  Mill'  simulation  of  nuclear  burst  phenomena.  This  requires  ;i 
detailed  knowledge  ol  the  rates  >u  mass,  momc'iiiiim,  and  energy  t.i.iiisfvi' 
between  the  ion- elect  roll  i In  id  and  the  neutral  fluid.  The  resistivity  of 
the'  medium  is  a ! so  recpiired  in  order  to  determine  the  decay  of  currents, 
perturbed  magnetic  fields,  and  st  ri.at  ions;  the  resistivity  is  directly 
related  to  the  rate  of  electron  momentum  loss.  Very  late- time  nuclear  burst 
simulation',  require  all  average*  collision  frequency  ol  ions  with  neutrals, 
and  this  collision  frequency  can  be  expressed  in  terms  of  the  rate  of  ion 
moment  urn  loss. 

line  major  approx  niiat  ion  that  we  have  made  in  ea  ieulat  in,;  the  trans- 
fer rates  due  to  collisions  is  that  we  have  considered  the  particle  velocity 
d i ' t r i bn  t tuns  of  cad.  species  to  he  Maxuel  i Hoi  t lanann  in  the  t I'ame  iiiO\  i in;  at 
the  mean  veloc  ity  ,,t'  that  species.  Ihts  approximation  is  a compromise:  it 

is  simple  enough  so  that  we  can  treat  many  cast's,  while  at  the'  same  time  it 
gives  sufficiently  accurate  results  for  use  in  high- alt  i t tide  simulations. 

Sect  ion  .!  derives  the  integrals  needed  to  cam  put  e the  collision;!) 
rales  of  momci'tum  transfer  (l.q.  J-  m>).  of  total  energy  transfer  (lq, 
and  of  internal  energy  transfer  (lq.  .!-  I0|  , and  the  rate  of  energy  loss  due  to 
inelastic  collisions  ilq.  11).  iliese  transfer  rales  are  expressed  in  terms 
of  tin-  women  t urn- coup  1 i ng  coefficient  s !lq.  .‘-.vil.  the  heat  - 1 rans  fer  coeft  i 
c i cut  li  (l.q.  and  the  inelastic'  energy  loss  integral  I.  (lq.  .’  xsi. 
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Wi  ion  > summai  i 7.cs  tin-  im  Miu'iH  Min  an-.!  int  ernal  energy  I ransfcr 
ratfs  needed  in  lim  -fluid  Mill*  siimilal  ion-.  siihso<|uciit  to  t !sc  first  Cow  set  - j, 

onds  aft  ■ r -i  inn  lour  I'urst  . Ilii  s 'avi  laii  i . sell  ionl. lined,  so  i e. tilers  who  i 

.no  }'l'l  n.  ill  Is  i ill  nvs  t til  in  results  i.m  skip  t In'  analysis  in  Sections  , a,  i 

.mil  I.  I iju.it  ions  ;.*JI  through  S show  tin-  simp  lost  version  oi  the  moment  inn-  4 

coupling  ami  lieu  t t run  a for  coefficient  - that  we  he  I I eve  to  he  idr.piut  . lor 

mu' I ear  hurst  s iimilut  i ons.  l.j.  o-.Mi  show-.  the  simplest  t'orin  of  the  electri-  1 

out  resist  i\  it}  uOeipiute  over  the  electron  t ompor.it  tire  range  of  0 . til  to 
In  eV.  lap..  V il  s|mn  l!u'  simplified  form  of  the  collision  I i etpienc  i es  for 
collisions  of  ions  Kith  neutrals,  of  electrons  tilth  iieuti.il'.,  and  of  elec- 
trons Mth  ions  t hut  .in'  needed  for  \er>  lute  time  s nmilut  toils . 

Purine,  the  initi.il  dclu  i s-a  i r interact  ion  when  the  particle  veloc- 
ities exceed  1(1°  iiiisec.  nni.nl  n collisions  are  more  important  l hull  elastic 
collisions.  Appendix  t summon  ~cs  i he  relevant  experimental  stopping  crons 
sect  tons  tluMi  i li  and  shots  s how  to  relate  rhese  to  the  moment  um- 1 runs  for 
cross  sis  t i on  ij  1.1'.,.  t Si  and  the  total  scat  tori  it;',  cross  section  o (laj, 
t-p  which  are  needed  to  compute  (lie  collisional  iiiomonttim  ami  energy  transfer 
rates  fiom  l.ijs.  ..  i.i  through  J II. 

Iw.i  fluid  Mill*  s nun  I at  iitits  are  cs.entiul  for  realistic  . i n.u  I at  i on  - 
of  lule-t  i me  mu  lour  hurt  plteiionu  na  . I hifoi  t im.il  i 1 1 . the}  arc  extreme!} 
complex  and  therefore  cannot  he  redone  easily.  Iluis,  i(  ; important  t ha  I 
the  momentum  and  energy  transfer  c oof fie i out s he  reasonahlx  accurate  te.g,, 
to  witliin  a ! 1 i.  iVe  heilCVe  t!ie  f.)  I I uw  i lie.  aiiuiys’s  altii  cross  SCI  t ion  dai.l  do 
provide  a reliahle  set  of  iocl  fuieitt  s,  I f the  reader  notes  any  errors  in 
this  uli.:l's|s  or  is  aware  ol  heller  cross  section  data,  we  would  appreciate 
hearing  from  him  so  iltat  improvements  can  he  made  as  soon  as  possible. 

Iii  help  the  reader  Keep  track  of  t ho  mat hemat  ica ! notation,  a 
bios sa ry  of  Symbols  is  given  in  Appendix  A. 
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SECTION  2 

THEORY  OF  MOMENTUM  AND  ENERGY  TRANSFER  BY  COLLISIONS 


2.1  SCATTERING  COLLISIONS 
(a)  Statement  of  the  Problem 


We  consider  two-bodv  (inelastic  or  elastic)  collisions  character- 
ized by  an  arbitrary  differential  scattering  cross  section.  The  colliding 
particles  may  be  atoms,  molecules,  ions,  or  electrons.  The  particles  of  the 
first  species  are  characterized  by  mass  m^ , uniform  number  density  Uj , 
mean  velocity  V^,  and  temperature  T^  (the  particle  velocity  distribution 
is  Maxwell- Boltzmann  in  the  frame  moving  at  velocity  V );  the  particles 
of  the  second  species  are  characterized  by  m7,  n^ , V,,  and  T?.  In  a 
collision  between  two  particles  of  different  species  moving  initially  at 
relative  speed  v,  each  particle  is  scattered  through  an  angle  G measured 
in  their  ccntor-of-mass  frame,  and  in  inelastic  collisions  the  pair  loses  an 
amount  of  kinetic  energy  c due  to  ionization,  excitation  followed  by  radia- 
tion, etc.  The  collisions  are  described  by  a differential  scattering  cross 
section  do(v,0)/d£2  in  the  centcr-of-mass  frame  (Reference  1).  The  inter- 
action potential  is  assumed  to  be  spherically  symmetric,  so  the  azimuthal 
angle  f is  not  needed. 

The  scattering  cross  section  and  the  inelastic  energy  loss  e 
depend  on  the  initial  and  final  states  of  the  particles,  ior  elastic  colli- 
sions (Section  3)  these  are  all  ground  states,  l-'or  inelastic  collisions 
(Appendix  C)  we  must  do  sums  over  states. 
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l.et  the  colliding  particles  have  incoming  velocities  v. , v and 
outgoing  velocities  v ^ , v,  tn  the  laboratory  frame.  (No  use  lower  case 
v's  to  denote  particle  velocities  and  capital  V's  to  denote  macroscopic 
velocities  of  gases.)  The  density  of  particles  of  the  first  species  in 
velocity  space  is  Maxwel 1- Boltzmann  in  the  frame  moving  at  velocity  V ^ ; 
wc  denote  these  Maxwel 1 -Boltzmann  distributions  by 


and  likewise  for  the  second  species.  Here  <S  is  the  Dirac  delta  function. 


Wo  take  the  rate  of  change  due  to  collisions  of  the  density  of 
momentum  of  the  first  gas  to  be 


~A—  = J "llvl  -V1Jtl(Vl)f2(V2)'VrV2l 


do(| Vj-v^l ,0j 


an 


dfld3v]d3v2 


(2-2) 


and  that  of  the  total  kinetic  energy  of  particles  of  the  first  species  to  be 

•'^V'i  ) f i r .2  2]  r ...  , v i r>‘  -y  ,ao^vi'v2^0J 

" I J mi[Vl  'V1  Jtl(vl)t2tV2,‘VrV2' 


At 


dQ 


dikf’v  ^d“'v7  , 


(2-5) 


where  pj 
Holt  zmaiiii 


= These  quantities  derive  from  the  collision  term  in  the 

-> 

transport  equation,  (dfj(Vj)/dt)  jj  for  example,  the  first  one 
times  this  collision  term,  integrated  over  all  velocities  v^. 


Our  terminology  concerning  the  forms  of  energy  of  a gas  follows  a 
scheme  used  by  some  workers  in  high-altitudc  phenomenology.  The  internal 
energy  of  a gas  is  the  translational  kinetic  energy  of  its  particles  in  the 
frame  moving  with  the  mean  velocity  V ^ plus  the  energy  of  rotational  ex- 
citation. (In  this  scheme,  the  ionization  energy  and  the  energies  of 
electronic  and  vibrational  excitation  must  be  treated  separately.)  This 
internal  energy  density  is  denoted  by  pi  The  "total  energy"  of  a gas 
is  the  sum  ot  its  internal  energy  plus  its  macroscopic  kinetic  energy  j p v^, 
We  assume  unless  otherwise  noted  that  changes  in  rotational  energy  are 
negligible,  and  so  we  refer  to  A / j p^2  \ simply  as  the  change  of  total 
energy  density  of  . ' first  gas.  ' 


\n  the  integrals  in  liqs.  2-2  and  2-3  the  final  velocity  v ' is 
•>  ->  * 
a function  of  v.,  v7,  e,  and  0.  We  will  change  the  independent  variables 

~r  ->■  _->  -> 

to  Vj,  v v - v.,,  c,  and  0 because  -he  differential  scattering  cross 
section  is  most  easily  expressed  in  of  these  variables. 


(b)  Change  of  Momentum  and  Kine J ergy  in  a Collision 


hirst  we  seek  to  o'- 
terms  of  the  (new)  indeoo 
these  species  by 


■ - final  particle  velocity  v^'  in 
les.  We  denote  the  reduced  mass  of 


U = iiyi^/  (iDj  +in2g 

the  ccntcr-of-mass  velocity  by 


(2-4) 


cm 


Vi  + "12v2 

nij  + m2 


(2-5) 


and  the  initial  relative  velocity  by 


V = V -V 
1 2 


(2-0) 


The  initial  velocities  are  related  by 


The  scattering  angle  0 in  the  center-of-mass  frame  is  defined 
by  the  diagram. 


The  azimuthal  angle  <i>  measures  the  angle  around  the  v axis.  We  assume 


that  the  differential  cross  section  is  independent  of  <|>,  so  we  can  average 


.>102 
Av^  and  T Cv^  ) over  <f>.  The  component  of  Av  perpendicular  to 

v therefore  averages  to  zero,  while  the  parallel  component  averages  to 


Av  7 = 


v • Av  -► 
v 


(v«v  -V  ) 


- - ^1  - — - ■ cosoj  v 

> 

or  by  the  conservation  relation  F.q.  2-9, 

/ 

V 

= - (l  - ./ 1 - COSo)  V 

^ \ ^ pv"  } 

-1 

'Hie  quantity  (Av)“  is  independent  of  (J> 

and  is  given 

(Av)^  = + v'-  - 2vv'cosO 

• v2(2  - 2L  . 2 

2e  fl\ 

y COS0 

UV”  / 
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(2-lla) 


(2-12) 


(2-13) 


If 


Vv£ 

■$k 


Using  those  in  relations  2-10  and  2-11  gives  the  two  desired 


averages, 


and 


<Av\>,  = - - n/1”  “2 cos0)  * 

ii W ■-  '"r*  (‘ - \[ 1 - cos0) 


(2-14) 


(2- IS) 


(c)  Evaluation  of  integrals 


Having  expressed  the  final  velocity  in  terms  of  the  initial 
velocities  and  0 and  e,  we  arc  ready  to  integrate  Hqs.  2-2  and  2-5 . 
First  wc  express  t he  product  of  the  velocity  distributions  in  terms 


of  the  \ iriablcs  and  v.  We  define  the  abbreviations 


, ml  Vm2  (VV 

tl\Vl;  “ (m1"+»2) 


-►  •>  -> 


v)  m(t2-t1)(v1-v2-V) 

(n-^Ti) 


W = « 


nj^tVj-Vj) 


- -V  /->  ->  - V->+»2’r,\ 

f2(v)  5 f2^v;  , 


(m 


T1T2  ) 
1 1 2+m2 1 1 } ' 


if  T > 0 or  T,  > 0, 


if  Tj  = 0 and  T_,  = 0, 
(2-16) 

(2-17) 


f7iv)/n7  is  the  normalized  distribution  of  relative  velocities  of  pairs  of 
particles  of  different  species.  Some  algebra  then  shows  that 


f.(v  V s})f2(v  'V  Sj-)-  V?i>f2(?>. 


(2-13) 


and  the  Jacobian  of  the  transformation  has  an  absolute  value 


ac’vj.v) 


^CVj.v,) 


(2-19) 


so  that  d'^v^d^v^  = d'\^d'\r. 


In  treating  t'  (v)  it  is  convenient  to  denote  the  macroscopic 
velocity  of  the  first  gas  with  respect  to  the  second  by 


V E V - V 
1 2 


(2-20) 

and  the  temperature  characterizing  the  distribution  of  relative  velocities  by 

(2-21) 


(m.T7+m,T  ) 
T = ” 


in  ~ (nij+m,) 


The  subscript  m denotes  a weighted  mean.  We  also  define  the  total  scat 
tering  cross  section  (Reference  1) 


Ot  (v)  .f  S 


dS 
dii 


(2-22) 


4ir 


and  the  momentum- transfer  cross  section 
q(v) 


/ (‘  ' f'  - J cos 6) 


(2-25) 
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This  cross  section  is  called  Q^c<1  by  McDaniel  and  Mason  (Reference  2, 
l>  158).  for  elastic  collisions  (f.  = 0)  we  denote  this  cross  section  by 
Mp(v) ; it  is  also  known  as  the  diffusion  cross  section. 


J 


liquations  2-2  and  2-3  then  become 


MpjVJ 


At 


— = -feliv\ye,(v)VqMuv\i\1d-\  , 


(2-24) 


% 


-vf  ) 


At 


* J'fl(v1)f2(v)v  v2-v*v1|q(v)  - ^ £Ot(v)  d\1d'\  . 


(2-25) 


The  integrals  over  d are  done  next,  and  afto''  some  algebra  on  the 
energy  equation  we  have 


A(PjVj) 


f2(v)q(v)vvd'>v  , 


f (\  - M CT^-T  ) 

|v  j L“  ~vsf-  • v + C^T^)  (v~'v‘v) 

- "l  ,«7  CP 2(v>)0tl 


The  last  integral  in  l.q.  1-11,  J f,  (vjo^  (v)  vd‘\ , is  the  moan 
frequency  of  scattering  collisions  of  a particle  of  the  first  species.  We 
will  use  the  closely  related  rate  coefficient  for  scattering  collisions  (of 
inelastic  energy  loss  c), 


Mvv> 


(v)  vdJv 


(2-28) 


We  need  two  other  integrals,  for  which  wo  use  the-  closely  related  rate 
coefficients  s and  h: 


s(Tn,V)  V ^ ~~  H ft'2(v)q(v)vvd‘*v 
hCT»*V)  /*/v)q(v)v(v2-v-V)dA 


(2-29) 


v . (2- SO) 


Note  t li.it  k and  h are  in  cm" /sec , while  s in  in  gm-cm''  / see.  We  will 
show  in  Section  2. Id  that  s can  he  interpreted  as  the  rate  coefficient  for 
momentum  transfer  or  "coupling  coefficient"  and  h as  the  rate  o efficient, 
for  heat  transfer. 

We  do  the  angular  integrals  using  spherical  coordinates  about  the 
V axis.  Integrating  first  over  the  azimuthal  angle,  we  note  that  in  the 
integral  of  liq.  2-29,  the  components  perpendicular  to  V integrate  to  zero. 

This  result  occurs  because  the  particle  velocity  distributions  of  each  species 
are  considered  to  be  Maxwell -Boltzmann  in  some  <' came . Wc  then  integrate  over  the 
cosine  of  the  zenithal  angle  and  change  the  iv  .lining  integration  variable  to 


\j  2k  (mjT-^in-Jj) 


For  brevity  wc  eliminate  the  mean  temperature  T in  favor 

1 in 

thermal  speed  parameter 


(2-31) 


tV;  related 


(2-52) 
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The  results  for  the  three  rate  coefficients  are 


- Y-f5  .if2)  e-«>V/U>‘ 


if  T > 0 aiul  V > 0 


A s(T  ,0) 
v in’  J 


if  V = 0 , 


^ vf,iq(v)+v  ii(v  i) 


if  T = 0 


'-(f-v/ur  -«>v/ur 

e - 0 


i f T > 0 and  V 
m 


<VVV)  = ' 


> 

l)  fro  (ligcV^  d(, 
Jr.  J ^ 


r v = o , 


Vrc>t  (V) 


if  i()l  “ o . 


(2-35)  * 


Those  rate  coefficients  art  continuous  functions  of  "f  (or  If)  ami  V.  The 

cross  sections  q needed  to  evaluate  the  integrals  (2-33)  in  (2-S!>)  are 

shown  in  figs,  l-l  to  1-8.  The  quantity  ca  can  be  related  to  the  cross 

) '• 

sect  ion  q by  co(  - v"q^u  j in ^m-,/ (m ( *111.,)  (f.q.  C-?  of  Appendix  C). 

(d)  General  Results 

The  time  rate  of  change  due  to  collisions  of  the  density  of 
momentum  of  the  first  gas  is 


ACPiV 


11.11  ,s i T . V I (\  -V 

1 .:  m 1 


(2-36).  . * 


and  that  of  the  total  energy  of  the  first  gas  is 


4?  pivi2) 


* ► > 

(.III  ( V j till  ,V  , ) • ( \ j-V, ) 

( III  I t 111  ( ) 


S (i  ,V) 
III 


+ - Tw^cksO>V)  . 


(2-37). 
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The  total  energy  of  each  gas  can  always  he  split  into  its  macroscopic 
kinetic  energy  aiul  its  internal  energy: 


At 

The  rate  of 

change 

of  the  i 

due  to  coll 

i s i ons  i 

s 

■'(§ 

-.v) 

At 

so  wo  have 

for  the 

rate  of 

collisions 

V 

At 


At 


[2- Mi) 


At 


l2-ol>) 


At 


1 


— s (T  , V ) V " «■  It  (T  , V)k  (T  -T  ) 
(ii^uuj  in  ' in  2 1 


ill. 


(liij+m.,)  L^s^1m’V^ 


]• 


(2-40) * * • 


The  first,  term  on  the  right  side  of  hq.  2-10  is  the  part  of  the 
macroscopic  kinetic  energy  of  the  gases  that  is  converted  to  internal  energy 
of  the  t’irsi  gas;  this  is  the  rate  of  internal  or  resistive  heating.  The 
second  term  is  the  heat  transfer  between  the  gases,  and  the  third  term  is 
the  energy  diverted  to  excitation  and  ionization.  We  refer  to  h as  the 
rate  coefficient  for  heat  transfer  and  to  s as  the  rate  coefficient  for 
momentum  transfer  or  simply  as  the  coupling  coefficient. 


Rates  for  the  second  gas  can  be  found  from  symmetry  by  inter- 
changing subscripts  1 and  2.  The  total  power  per  unit  volume  transferred 
by  collisions  from  the  total  energy  of  both  gases  to  other  forms  of  energy  is 


At 


At 


11,11  ,ck  (T  . V ) 
12  s'  in 


(2-4  1 ) 
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I or  tin.'  spec  i a I i-asr  of  elastic  collision:;  hot  noon  two  gases  having 
no  iv  I at  l vr  volooity  t<:  = 0,  \ - if),  l.q,  agrees  with  hvsloge’s  result 

(Reference  a,  l.i|.  Ji'I. 

i.Z  CHARGE- EXCHANGE  REACTIONS 
(<i)  Introduction 

While  elastic  seat  ter  ini',  is  the  most  important  means  of  trans- 
ferring momentum  and  energy  between  ions  and  neutrals  at  low  speeds,  charge 
exchange  is  also  important  at  relative  pavtiele  speeds  tihove  in'*  era/ see.  In 
this  report  we  concentrate  on  symmetric  charge  exchange,  in  which  the  products 
are  the  same  as  the  reactants: 


\ * X • X • X 


U'-4.)> 


It'  each  reactant  and  product  is  in  its  ground  state,  this  reaction  has 
"ero  energy  defect  and  is  called  resonant  (or  elastic)  charge  exchange. 

Such  resonant  reactions  can  have  unusually  large  cross  sections  (see 
figures  i-l  to  4-4)  and  are  among  the  most  important  reactions  at  times 
after  a few  seconds  after  a high-altitude  nuclear  burst.  Symmetric  charge 
exchange  is  much  like  a scattering  process,  and  we  will  treat  it  hy  a minor 
modification  of  the  scattering  theory  of  Section  J.l. 

Asymmetric  charge  transfer  is  a chemical  reaction  in  our  view 
and  is  treated  differently.  A treatment  comparable  to  Section  . I is 
given  in  Appendix  1!,  subsections  (al  to  (d),  and  a simplified  treatment 
is  given  in  subsection  (e),  lixperi  mental  cross  sections  for  some  asym- 
metric charge -exchange  reactions  are  plotted  in  Section  4 for  comparison 
with  the  others.  We  do  not  want  to  become  enmeshed  here  in  choosing  which 
chemical  reactions  are  important,  though,  so  wo  have  not  evaluated  rate 
coefficients  for  asymmetric  charge  exchange  in  Section  5. 


(b)  Symmetric  Charge  Exchange 


Symmetric  charge  exchange  can  he  treated  within  the  iramewo rk 
of  Section  2.1,  he  simplify  this  slightly  by  assuming  that  such  reactions 
are  forward  scattering  mui  have  ivj  inelastic  loss  or  gain  of  kinetic  energy 
(that-  is,  dn/dw’  is  neg  i ig il>le  unless  0 « i and  c = 0], 


Symmetric  charge  exchange  can  he  considered  an  ordinary  scatter- 
ing collision  by  relabel  inr,  the  part  icles  after  the  encounter.  An  encounter 
that  is  forward  scatter  inr,  in  the  cvnter-of-itiass  frame  is  backseat  t «»>*»  t»j* 
m this  alternate  description.  The  effective  moment  urn- 1 vans  for  cross 
section  in  this  alternate  description,  denoted  by  an  asterisk,  is  given  by 
kgs . 2-23  and  2-22  »*ith  u , « and  t » tt: 


g*  iv  i « 2ox(v) 


(2-43) 


This  point  is  elahoratod  in  Reference  , Section  5-7-A.  divou  o (v)  one 
can  use  this  relation  in  Lqs.  2-33  and  2-34  to  get  a coupling  coefficient 
s*  ami  heat  - 1 runs  for  coefficient  It*  due  to  symmetric  charge  exchange; 
these  must  he  added  to  the  coupling  coefficient  s and  heat- 1 ransi'er 
coefficient  h due  to  scattering  collisions  (see  kg.  2-48). 

2.3  RELATION  10  HIE  EQUATIONS  Of  FLUID  DYNAMICS 
(a)  Momentum  Balance  Equations 


We  show  hero  how  our  results  for  elastic  scattering,  and  syiume 
charge  exchange  relate  to  the  momentum  and  energy  balance  equations  of 
fluid  treatment.  Inelastic  scattering  is  not  considered  here;  at  times 
l see  following  a nuclear  hurst  when  the  two-tluid  equations  are  a good 
mation,  inelastic  scattering  is  unimportant  except  for  the  energy  t vans 
tween  electrons  and  the  vibrational  excitation  of  X,  molecules. 

file  neutral  fluid,  indicated  by  subscript  n.  is  composed  of 
species  \°,  j ■■  1,2,...,  so  that 

n - 

n ? J 


t v i c 
a two* 
aft  el' 
approx i 
for  be- 


(2  4 1) 


jUJuljj* to  llkt,,, 


similarly  the  ions  are  composed  of  singly  charged  species  X^,  k = 1,2,..., 


so  that 


ni  = IIn0C)  • 

k K 


(2-45) 


We  assume  that  all  neutral  species  have  the  same  temperature  T and  gas 
velocity  and  all  ions  have  the  same  temperature  and  velocity 

V. . The  electrons  have  density  ng  = ir  , temperature  T , and  velocity 

(2-46) 


-v  -y 
V = V. 
e i 


cJ/en. 


for  cacli  pair  of  a neutral  species  j and  an  ion  species  k, 
Eq.  2-21  defines  a weighted  mean  temperature 


_ (mj'1Vi\Tn) 

jk  (in  j +l\ ) 


(2-47) 


If  either  species  is  an  electron  the  weighted  temperature  is  simply 


The  rate  of  momentum  transfer  due  to  ion-neutral  collisions  is 
a sum  of  terms  of  the  form  of  the  right  side  of  Eq.  2-36.  For  each  pair 
of  a neutral  species  j and  an  ion  species  k,  Eq.  2-33  defines  a coupling 
coefficient  sjg  H jg > I I ) due  to  scattering.  The  ion-neutral  coupling 
coefficient  s ^ is  given  by 

"„v„i 

+ En(X°)n(X*)s*((Tn+T.)/2,|Vn-Vii)  , (2-48) 

k 


where  the  s*  are  coupling  coefficients  uue  to  symmetric  charge  exchange 
of  ions  of  species  k with  their  parent  atoms  or  molecules  (Section  2.2b). 
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Similarly  the  electron-neutral  coupling  coefficient 


s 

ne 


i s 


11  s 
n ne 


£n(X°)s.  (T  ,|tf  -Vl) 
V ijee'n  e ‘ ' 
} 


(2-49) 


where  we  have  assumed  that  the  neutral  temperature  is  small  compared  to 
(m./m^JT^.  The  ion-electron  coupling  coefficient  is 


n.  s . 

l le 


£n(X!^Ske(Te*cJ/eni)  * 
k 


(2- SO) 


s.  and  s arc  closely  related  to  the  electrical  resistivity  (Kq . 5-22). 
le  ne 


The  momentum  balance  equations  for  the  neutral  fluid,  for  the  ion- 
electron  plasma,  and  for  the  electrons  can  be  written  as  (see  Glossary  of 
Symbols  in  Appendix  A) 


3(0  V ) 

n n 

3t 


+ V‘(o  • V ) a - VP  + p g - n n.s  . (\/  -V.)  - n n s (V  -V  ) 
n n n n n n l m n iJ  n e ne  n e 


(2-51) 


3 (Pi  V. ) 


IT-  * Hpxv  = - voyiy  + 3 x b ♦ (p.4.pe)g 


0 = -VP  - en 
e c 


+ n n. s . (V  -V. ) + n n s (V  -V  ) . 

n i ni  n l n e ne  n e * 

/ V xg\ 

/■►  el  ^ ->■  *► 

11;  + — - — / + p g + n n s (V  -V  ) + n.n  s.  (V.-V  ) 

\ cl  n e ne  n e l e iev  l c 


(2-52) 


(2-53) 


We  have  ignored  viscosity  and  electron  inertia.  One  can  eliminate  V in 
terms  ot  V\  and  J by  liq.  2-46. 


(b)  Energy  Balance  Equations 


for  the  energy  balance  equations,  Uq.  2-34  defines  a heat-transfer 
coefficient  h^  (1\^ , IV^-V^  | ) for  each  pair  of  a neutral  species  j and 
an  ion  species  k.  The  ion-neutral,  electron-neutral,  and  ion-electron  heat- 
transfer  coefficients  h h , and  h.  are  defined  as  in  liqs.  2-48  to  2-50. 

XI 1 II O 1C 
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The  energy  balance  equations  also  include  frictional  or  resistive 

heating,  which  is  the  sum  over  all  pa:.  of  species  of  the  first  term  on 

the  right  side  of  Eq.  2-40.  For  each  pair  of  species  the  rate  at  which 

collisions  between  species  1 and  2 convert  their  total  macroscopic  kinetic 

2 

energy  into  internal  energy  is  n^n2s(lm,V)V  , and  each  species  receives  a 
fraction  of  this  that  is  inversely  proportional  to  its  particle  mass.  The 
fraction  d 


(n) 


ru 


of  the  ion-neutral  frictional  heating  that  goes  into  internal 


energy  of  neutrals  is  given  by 

n n.s  .d^  = X.  n(X?)n (X^ ) 

n l m ni  • V j k' 
1 k 


\ 


-v.i) 


(nw+i«k)  jkL  jk'1  n i1 
£ »(X°)n(X*)  i s*(i  T„  . i Ti.|?n41|)  , (2.53a) 


recall  that  j and  k specify  neutral  and  ion  species,  respectively.  The 

f _ •» 

corresponding  fraction  that  goes  to  ions  is  1 - d\  . For  atmospheric  ions 

and  neutrals  considered  in  Sections  3 to  5 this  fraction  is  in  the  range 

0.3  < d*'.  < 0.7.  We  also  define  a "fractional  coupling  coefficient"  for 
ni 

electron  neutral  collisions 


d ^ by 
ne  1 


n s d(n)  = 23n(X?) 
n nc  ne  ~ \ 


7 - — r S . (T  , 

(m . -Hn  ) l e e 
j e 


V -V  j) 
1 n o 


(2-54) 


and  similarly  for  d) 
tlie  electr 
and  2-65) . 


(i) 


The  fractions  d 


(n) 


(i) 


and  d : 

1C 


arc  of  the  order  of 


ic  nc 

the  electron  mass  divided  by  the  mean  mass  of  neutrals  or  ions  (see  Hqs.  2-63 


With  tiie  collisional  terms  expressed  by  means  of  Cq.  2-40  with 
e = 0,  the  energy  balance  equations  for  the  neutral  fluid,  for  the  ions, 
and  for  the  electrons  are 

‘HP,,!,.)  -+■  + (nl ->■  -*■  2 

2-2-  + V • (p  I V ) = - P V • V + n n.s  .dr  /(V  -V.) 

3t  n n n n n n l ni  ni  n i' 

+ n n s d^n)  (V  -V  )2  + n n.h  .k(T.-T  ) + n n h k (T  -T  ) , 

n c ne  ne  n eJ  n i ni  l lr  n e ne  e n 


26 


(2-55) 


i-i-  + V • (p.  I.  V.  ) = -P.7-V.  + n n.s  , (l-d^HV  -V.) 

dt  ^1  1 i i n l nr  ni  J K n i 

+ n.n  s.  d.(l)  (tf  -V  )2  - n n.h  .k(T.-T  ) - n.n  h.  k(T.-T  ) , 

l e le  re  i e n l ni  i n i c te  v i e 


(2-56) 


3(n  4 kT  ) 

c 1 e 

at 


(nc  I " 


l‘  V • V + n n s (V  -V  )' 
e e n e ne  n e 


n.n  s.  (V.-V  y - n n h k(T  -T  ) + n.n  h.  k (T. -T  ) . 

l e iev  i e'  n e nc  e n i e le  l e' 

(2-57) 


We  have  ignored  thermal  conductivity  and  viscosity,  which  are  not  treated  in 
this  paper  although  they  are  important  processes  above  300  km  altitude.  In 
liq.  2-57  we  have  assumed  that  1 - d^1^  - 1 and  1 - d.^  « 1 . In  the  approxi- 
mation  that  \T  » V , the  second  and  third  terms  on  the  right  side  of  this 
equation  are  the  rate  of  resistive  heating  qJ2,  where  n is  the  electrical 
resistivity  (liq.  5-22).  For  nuclear  burst  simulations  that  include  an 
vibrational  temperature  T , liq.  2-57  should  also  include  heat  transfer  pro- 
portional to  (T  -T  ) . 
i v e 


A special  case  provides  several  simplifications  to  the  rate 
coefficients.  Our  numerical  results  in  Section  5 are  for  the  special  case 
of  low  relative  velocities: 
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2k  T 


in. 
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2kT . 
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(2-58) 


|V  -V  I2  « 2kT  /m  . and  similarly  for  the  ion-electron  velocity.  From 
1 n e 1 e e 

liq.  2- 34b  the  heat-transfer  coefficients  become 


n n.h  . 
n l in 
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j k J 
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(m.+m.  ) 
j k 


£n(X°)n(X*) 
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(2-59) 
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We  - £ • O’*®’ 

• (2"61) 

k 

Electron-neutral  scattering  is  important  mainly  at  altitudes  below  150  km, 
where  the  neutrals  are  primarily  N?.  From  Eqs.  2-60,  2-49,  and  2-54  one 
finds 


"no  * 3 ' (2-62> 

ill /iii(N2)  = 2.0  * 10-5  . (2-63) 

Similarly  if  most  ions  have  a mass  near  that  of  0*,  one  finds  from  Eqs.  2-61 
and  2-50 

h.  ^ 3s,  /m(0+)  , (2-64) 

le  le  v 

d^)^me/m(0+)  = 3.4  x io'5  . (2-65) 

Under  the  above  assumptions  the  ion-neutral  fractional  coupling  coefficient 
defined  in  Eq.  2-53a  can  be  approximated  as 

d(n^  0.4  , (2-66) 

ra 

accurate  to  ± 20  percent. 


SECTION  3 

SCATTERING  CROSS  SECTIONS  FOR  ATMOSPHERIC  SPECIES 


To  use  the  preceding  theory  in  two-fluid  simulations  of  high- 
altitude  nuclear  explosions,  one  needs  momentum-transfer  cross  sections  for 
scattering  collisions  as  a function  of  relative  speed  or  of  energy.  One 
needs  to  consider  ion-neutral,  electron-neutral,  and  electron-ion  scattering 
for  principal  atmospheric  species.  Also  one  needs  charge-exchange  cross 
sections,  which  will  be  presented  in  Section  4. 

IVe  consider  the  species  N,  0,  N , 0 , NO,  lie,  N , 0+,  N +,  0,+  , 

+ + **  Z Z 

NO  , He  , and  electrons.  We  concentrate  on  elastic  collisions  (e  = 0)  in 
the  energy  range  from  0.005  eV  to  roughly  30  eV  for  electrons  or  tens  of 
keV  for  ion-neutral  collisions.  At  higher  energies,  inelastic  collisions 
dominate;  ion-neutral  inelastic  collisions  are  considered  in  Appendix  C. 

We  consider  electron-neutral  scattering  using  experimental  data,  electron- 
ion  scattering  using  theory,  and  ion-neutral  scattering  using  mainly  theory 
or  semi-empirical  methods. 

3.1  ELECTRON-NEUTRAL  SCATTERING  BELOW  100  eV 

For  scattering  of  electrons  from  N,  0,  N?,  0?,  and  NO,  we  have 
used  the  momentum- transfer  cross  sections  recommended  by  Phelps  (Reference  4). 
For  the  molecules  his  ci'oss  sections  extend  over  electron  energies  from 
0.005  to  100  eV.  He  tabulates  the  frequency  of  momentum- transfer  collisions 

3 

per  unit  density  of  molecules  (cm  /sec);  we  have  divided  by  the  electron 
speed  to  obtain  a cross  section.  Figs.  4-3  to  4-4  show  these  cross  sections. 
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The  ci'oss  sections  come  from  various  types  of  experiments.  The 

e + N?  cross  section  comes  from  an  analysis  of  dc  swarm  experiments,  and 
* * 
Phelps  estimates  ±20  percent  accuracy  or  better.  The  e + NO  cross  section 

is  based  on  an  analysis  of  dc  mobility  data  and,  for  energies  below  0.02  cV, 
on  theory  for  a polar  molecule;  the  accuracy  is  a factor  two.  The  e + 0^ 
cross  section  below  0.1  eV  is  from  microwave  measurements  and  is  considered 
accurate  only  to  order  of  magnitude;  at  higher  energies  it  is  from  beam 
and  swan  experiments  and  is  good  generally  to  ±20  percent.  The  e + N and 
e + 0 cross  sections  are  order-of-magnitude  estimates  based  on  total  scatter- 
ing cross  sections  from  electron  beam  experiments;  these  cover  a more  re- 
stricted range  of  electron  energy. 

For  e t lie  scattering  we  have  used  the  momentum- transfer  cross 
section  recommended  by  Itikawa  (.Reference  5),  which  covers  election  energies 
from  0.01  to  10  eV.  This  is  based  on  dc  swarm  experiments  and  should  be 
accurate  to  better  than  20  percent  (Reference  6).  This  is  plotted  in  Pig.  4-1. 

3.2  ELECTRON- ION  ELASTIC  SCATTERING 


We  consider  electron-ion  collisions  at  low  energies,  where  in- 
elastic processes  are  not  yet  important.  For  electron  energies  up  to  10  eV, 
we  need  consider  only  elastic  (coulomb)  scattering.  For  a coulomb  force 
between  a singly- charged  ion  and  an  electron  we  use  the  Rutherford  differ- 
ential scattering  cross  section  (Reference  1) 


do,  Q,  e 
dfl(V,6)  “2.2 


4u“v  sin  (0/2) 


(3-1) 


For  e + NO  at  20  eV  energy  we  use  2.95 
of  Phelps'  tabulated  value. 


x 10 


-7 


cnr/sec  instead 
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The  reduced  mass  is  1.1  * m , but.  we  have  not  used  this  in  order  to  retain 

e 

the  strong  advantage  of  having  the  results  symmetric  between  ions  and 
electrons,  in  the  Debye  approximation  the  force  is  assumed  zero  beyond 
one  Debye  shielding  length  A^  because  the  ion  charge  is  shielded  by  free 
electrons . 

The  minimum  centcr-of-mass  scattering  angle  is  then  (Reference  7, 
hqs.  8-34,  8-35) 


0 = maximum  / , 

m U^v2  Vv/ 


(3-2) 


which  gives  a finite  momentum- transfer  cross  section  (from  Eq.  2-23), 


. . 4ire  , 2 

Vv)  = ~2 ~4  *n  0~ 

u v 111 


(3-3) 


The  weak  v- dependence  of  £n(2/0  ) is  ignored  by  substituting  in  Eq.  3-2 
a typical  relative  speed. 


/3kT 3kT. 

/,  = /— -£■  ♦ - 

t yj  me  m. 


+ V 


(3-4) 


where  T and  T\  arc  the  electron  and  ion  temperatures  and  V is  the 
ion-electron  macroscopic  relative  speed,  given  by  liq.  3-10.  Kc  abbreviate 
the  coefficient  in  the  cross  section  by 


A - 


4ne 


(3-5) 


= 8.06  x io1'  (2/0  ) cm(l/sec4 


The  factor  2n(2/0  ) depends  weakly  on  ne,  T,  and  V and  is  plotted  for 
V = 0 in  Reference  7,  Figure  8-6,  In  Section  5 we  use  fcn(2/0  ) = 12. 


This  Debye  analysis  is  valid  only  when  a sphere  of  radius  equal 
to  the  Debye  shielding  length  contains  a statistically  large  number  of 
electrons  (Reference  7,  liq.  8-18). 


The  coupling  coefficient  s^(T^,V)  and  heat-transfer  coefficient 

h.  (T  .V)  due  to  elastic  scattering  of  electrons  on  ions  are  calculated  by 
iev  m'  ' 

integration  from  Eqs.  2-33  and  2-34.  The  integrands  involved  have  removable 
singularities  at  f,  = 0 and  can  be  integrated  by  parts  . The  results  can 
be  expressed  in  terms  of  the  error  function, 


erf(x) 


(2//F)  f o' t 


dt 


(3-6) 


The  results  are 
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(3-7) 
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2 „3 


(3-8) 


where  we  have  again  eliminated  the  mean  temperature  T in  favor  of  tlie 
related  thermal  speed  parameter  (L-q.  2-52) 


2kT  2kT. 

U = + 

V ,ne  "i 


(3-9) 
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Usually  the  ion 
U * (2kTe/mc)\ 


temperature  is  small  compared  to  (•1Vi/me)'l'e*  so 
The  ion-electron  relative  velocity  is  x'eloted  to  the 


current  density  by 


•f 

V 


= cJ/ene 


(•■5-10) 


The  rates  of  change  due  to  ion-electron  collisions  of  the  momentum 
and  internal  energy  densities  of  the  ions  are  given  by  Hqs.  2-36  and  2-40: 


ACPiV 

At 


n.n  (V.-V  ) s.  (T  , I V.-V  |)  , 

i e v i e ie  in  1 j.  e 


(3-11) 


At 


tin  _ ->  . 

7 — 7-  IV. -V  |-  S.  (T  ,|v.-V  I) 

(nu+m^)  | l e|  ie  m 1 l e1' 

+ k(T  -T. )h.  (T  , I V.-V  I) 
v e te  m l i 


(3-12) 


The  rates  for  electrons  are  found  by  interchanging  subscripts  i and  c. 


In  the  static  case  ( |v. -V^  j=D) , the  second  term  on  the  right, 
side  of  li<j . 3-12  is  the  rate  of  ion-electron  heat  transfer;  when  liq.  3-8 
with  3-5  and  3-9  is  substituted  this  term  agrees  with  the  long-known  form 
(Reference  7,  F.q.  9-60;  Reference  8,  Fqs.  5-30,  5-31). 


The  first  term  on  the  right  side  of  Fq.  3-12  is  also  familiar. 
In  the  relation  for  electron  temperature  change,  which  is  found  from  liq . 
3-12  by  interchanging  subscripts  i and  e,  this  term  can  be  considered 
in  the  c.ise  of  a fully  ionized  plasma  with  equal  temperatures  as  the 
electron-ion  resistive  heating  qj-.  If  V <<  U so  that  liq.  3-7b  can  be 
substituted,  this  electron-ion  resistivity  agrees  with  the  resistivity 


q . in  a woak  electric  field  and  a strong  perpendicular  magnetic  field, 
given  by  Spitzer  (Reference  8 , liq . 5-42)  and  l.ongmire  (Ueforence  7,  liq. 
U)-63) . 


The  electron-ion  coupling  coefficient  of  l-q.  3-7  is  essentially 
independent  of  ion  mass.  It  can  be  identified  with  t ho  above  electron- 
ion  resistivity  by 


.0)  " (c/O-  M,  • 


(3-13) 


-3/-> 

and  has  a temperature  dependence  of  T . 


for  relative  ion-electron 

1.  i ...  ^ T O / — 


speeds  V « U the  heat-transfer  coefficient,  h.  varies  as  ra.  T ' , 

to  l e 

At  higher  speeds  both  s,  and  h.  decrease;  the  function  in  Table  3-1 
* 1 10  10 

is  useful  here. 


Wc  have  checked  one  other  special  case,  for  a beam  of  ions  imping- 
ing on  stationary  electrons  (i\  = 0 and  = 0),  hq.  2-38  with  substitutions 

2-39,  3-11,  3-12,  3 -7a,  and  3-S  agrees  with  Reference  7,  liq.  9-55. 

-x2 

Table  3-1.  The  function  f(x)  = erf(x)  - (2//fr)xe  * . 


X 

smal  1 

‘ 0.2  

0.4 

0.6 

0.8 

1.0 

f(x) 

(4/3/n t)x3 

0.00586 

0.0438 

0.132 

0.266 

0.428 

X 

1.2 

1.4 

1,6 

1.8 

2.0 

2.5 

f(x) 

0.590 

0.730 

0.837 

0.910 

0.954 

0.994 

3.3  ION-NtUTRAL  AND  NEUTRAL-NEUTRAL  ELASTIC  SCATTERING 

(a)  Low-Energy  Limit  — the  Polarization  force 

We  consider  elastic  collisions  of  singly  charged  ions  with  atoms 

or  nonpolar  molecules  at  low  energies,  below  about  0.1  eV  in  the  cont.er-of- 

mass  frame.  The  dominant  long-range  force  between  these  particles  is  the 

attractive  force  between  the  ionic  charge  +e  (or  -c)  and  the  electric  dipole 

moment  induced  by  it  in  the  molecule.  As  long  as  the  separation  r between 
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the  ion  and  molecule  exceeds  5 * 10  cm,  their  potential  energy  is  that 
due  to  this  polarization  force, 

Hr)  = 


*> 

o~« 

2r4 


(3-14) 


where  u is  the  polarizability  of  the  molecule. 


We  use  a classical  theory  of  this  interaction,  for  the  errors 
caused  by  a classical  treatment  seem  to  be  mild  (Reference  9).  The 
momentum- transfer  cross  section  is  given  by  (Reference  2,  Eg.  5-7-48; 
Reference  11,  pp.  457-458) 


q(.(v)  = 2.210  nc 


MW  i 


m.m, 
i k 


(3-15) 


This  is  the  polarizability  limit  of  the  classical  l.angevin  theory  (Reference 
12).  Table  5-2  gives  the  polarizabilities  of  some  atmospheric  gases,  and 
figures  4-1  to  4-8  show  the  related  cross  sections. 


The  coupling  coefficient  sjg(-i'm>v)  and  heat -transfer  coefficient 

h.,  ff  ,V)  for  elastic  scattering  of  ions  on  molecules  bv  the  polarization 
,lk  m . . 

force  are  known  to  be  independent  of  temperature  when  the  gases  are  static 


Table  3-2.  Angle-averaged  dipole  polarizability  and  electric  dipole  and 
quadrupok  moments  of  some  atmospheric  atoms  and  molecules  in 
their  ground  states  (Reference  2,  Appendix  II). 


Polarizability 
Neutral  a 

Molecule  j 

i 1 0"  *" 4 an 


0.205 
1.13 
0.77 
1.76 
1.60 
1 .70 


Di pole 
moment, 

Kf18  statcoul-cm 

Quadriipole 

moment, 

10  c statcoul-cnr 

0 

0 

0 

0 

0 

0 

0 

-1.52 

1 0 
1 

-0.39 

| 0.153 

-1.8 

(V  • 0);  this  result  remains  true  when  the  gases  are  moving.  From  liqs.  2-33 
and  2-34  the  rate  coefficients  become 


s..  (T  ,V)  1 2.21 0 rc 
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(3-16) 


(3-17) 


Table  3-3  lists  values  of  these  rate  coefficients  for  various  ion-molecule 
pa  i rs . 

The  NO  molecule  lias  a small  intrinsic  electric  dipole  moment  (see 
Table  3-2).  This  dipole  moment  is  essential  to  e + NO  scattering  at  low 
energy,  as  noted  in  Section  3.1.  For  scattering  of  an  atomic  ion  from  NO, 
this  dipole  moment  presumably  causes  a long-range  force  whose  potential  varies 
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3-3.  Couplinq  coefficient  and  heat- transfer  coefficient  due  to 

polarization  force  for  various  atmospheric  ion-neutral  pairs. 
Those  are  calculated  from  Equations  3-16  and  3-17. 
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Less  accurate,  because  of  the  intrinsic  <1  i polo  moment  of  NO.  See  text. 
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with  separation  like  r ; this  potential-energy  term  also  depends  on  the 
orientation  of  the  NO  molecule,  for  the  polar  gas  CO,  however,  measure- 
ments of  the  mobility  of  alkali  ions  in  CO  unexpectedly  agree  with  the 
prediction  based  on  the  polarization  force  (Reference  13;  Reference  10, 

Section  9-9-F).  This  agreement  may  be  because  the  dipole  moment  of  CO 

“18  _ 2 

is  very  small  (0.11  * 10  statcoul-cm)  or  because  the  r potential 
energy  vanishes  when  averaged  over  all  orientations  (Reference  10,  p.446). 
We  will  therefore  use  the  polarization-force  cross  section  and  rate  coeffi- 
cients for  NO,  as  shown  in  Figure  4-8  and  Table  3-3.  We  similarly  ignore 
the  electric-quadrupole  interaction  potential  (which  varies  as  r 3) 
between  an  ion  and  an  N?  or  0^  molecule. 

(b)  Scattering  Below  20  eV 

We  consider  elastic  atom-atom  or  ion-molecule  collisions  for 
center-of-mass  kinetic  energies  below  20  eV.  For  atmospheric  species  this 
implies  a relative  speed  less  than  2 x 106  cm/sec  and  a distance  of  closest 

_ g 

approach  of  about  1 x 10  cm  or  more. 


For  two  atoms  or  nonpolar  molecules,  the  interaction  potential 

energy  4>(r)  has  an  attractive  r"6  dependence  at  long  range  (the  London 

dispersion  or  induced-dipole- induced-dipole  interaction)  and  a repulsive, 

approximately  exponential  dependence  on  r at  short  range.  We  use  the 

-12 

Lennard-Jones  (12-6)  potential,  the  r term  being  a convenient  power-law 

approximation  to  the  short-range  repulsion.  For  an  ion  and  an  atom  or  non- 

-4 

polar  molecule,  the  dominant  long-range  interaction  is  the  attractive  r 
polarization  potential  considered  earlier;  for  this  case  we  use  a (12-4)  or 
(12-6-4)  potential.  We  express  the  potential  energy  in  terms  of  the  depth 
£ of  its  minimum,  the  separation  r at  which  tliis  minimum  energy  -e 
occurs,  and  a dimensionless  parameter  y measuring  the  relative  strength 
of  the  r 6 and  r 4 energies: 


*(r) 


(i 


12 

. . 1 rm  \ 
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(3-1 3) 


The  case  of  Y = 1 is  the  Lennard-Jones  (12-6)  potential,  and  that  of 
Y = 0 is  the  (12-4)  potential. 

Table  3-4  shows  the  assumed  form  and  the  parameters  of  the 
potential  for  several  pairs  of  atmospheric  species.  We  have  excluded  ions 
in  their  parent  gases  because,  at  room  temperature  and  above,  charge 
exchange  dominates;  a classical  calculation  would  be  dubious  and  the  re- 
sulting elastic  momentum  transfer  negligible  at  energies  above  0.03  eV. 
for  ion-neuti'al  interactions  we  choose  one  parameter  by  requiring  that, 
the  r~4  term  of  l-q.  3-18  match  the  polarization  potential,  Rq.  3-14; 
for  example,  for  N+  + 0 we  use  the  polarizability  a of  0 to  calcu- 
late Y by 


Y = 1 - 


(3-19) 


For  the  bottom  three  pairs  in  Table  3-4,  Mason  (Reference  14) 
has  estimated  ion  mobilities  for  the  temperature  range  0 < kT  < 0.22  eV. 
The  mobility  (cm^/statvolt- sec)  at  standard  density  of  a species  of 

singly-charged  ions  in  a weak  electric  field  is  related  to  our  coupling 
coefficient  s.,  for  static  gases  (V  * 0)  by  (Reference  12,  Section  5) 


= jk(r,0)  + s* (T,0)  - n K 
J o o 


(3-20) 
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where  n = 2.69  x 10  molecuics/cm  is  the  density  of  neutrals  at  standard 
o 

conditions  (0°C  and  i atm).  For  the  last  three  systems  in  Table  3-4  Mason 
used  the  mobility  of  the  H+  + H„  or  N0+  + N_  system  as  a model  and 
scaled  it.  An  equivalent  procedure,  which  we  adopted,  is  to  assume  that 
the  potential  energy  of  the  system  is  given  by  Kq.  3-18  with  £ and  Y 
being  the  same  as  for  Mason’s  model  system  and  r being  determined  from 
liq.  3-19  and  the  polarizability  a of  the  actual  system,  'liras,  for  these 


*7tt?** *~-— ‘ I 1 1 ~n  r i I I I ■ ■■■  I I,,-  . 


Table  3-4. 


Characteristics  of  potential  energy  functions  of  pairs  of 
atmospheric  species  for  separations  greater  than  1x10“®  cm. 


Collidi  ng 
pai  r 

Nature  of 
interaction3 

Form 
of  4>  ( r ) 

e 

eV 

rm 

10'8  cm 

y 

N + N 

val . 

12-6 

9.84c 

O 

o 

N + O 

val . 

12-6 

6.64c 

1 . i 5C 

0+0 

val . 

12-6 

5.22c 

1 . 22c 

He  + He 

rep. 

12-6 

0.00088f 

2.9f 

n2  + N? 

rep. 

12-6 

0.00749 

4. 29 

V °2 

rep. 

12-6 

0.00869 

3.99 

N+  + 0 

val . 

12-6-4 

4. 57c 

1 . 23c 

0. 65e 

0*  + N 

val . 

12-6-4 

10. 97c 

1.06c 

0. 61e 

N+  + N2 

val . b 

12-6-4 

~2.7b’d 

1.83e 

~0.72b,d 

0+  + 02 

val.b 

12-6-4 

~2.7b,d 

1.79e 

~0.72b’d 

02  + 0 

rep. b 

12-4(?)b 

~0.11b 
0. 32c 

■ . 

2.41e 

rep.  = short-range  repulsive  force 

val . = short-range  attractive  valence  force 

Mason,  1970  (Reference  14) 

Gilmore,  1972  (Reference  15) 

Mason  and  Vanderslice,  1959  (Reference  16) 


m 


i 


1 

■J 

if 


" -‘rl 

m 


e Calculated  from  Eq.  3-19  by  using  the  polarizability  of  the  neutral 
species  in  Table  3-2  together  with  e and  r (ory)  from  this  table. 
For  the  (12-4)  potential,  y is  zero. 

f Hirschfelder , et  al.,  1964  (Ref.  20).  Kestin  and  leidenfrost  (Ref.  21) 
give  r = 2.48  * 10"®  cm  and  a larger  energy,  e = 0.0060  eV. 
n m 

9 Hilsenrath,  et  al , 1955  (Reference  22) 
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systems  our  coupling  coefficient  s.|_  (T^.V)  corresponds  exactly  to  Mason's 

mobility  but  extends  to  higher  temperature  and  to  nonzero  relative  gas  vel- 
. *>  -*■ 

ocity  V = 


lVvil 


for  0^  + 0,  Mason  made  an  educated  guess  that  the  short-range 

force  is  repulsive.  The  dissociation  energy  given  by  Gilmore  (Reference  15), 
0 + -*  0 ' + 0 - 0.32  eV,  supports  this  assumption.  Mason  also  estimated 

-f 

the  mobility  of  NO  + 0 based  on  a short-range  repulsive  force,  but  this 
assumption  appears  wrong  in  view  of  Gilmore's  dissociation  energy. 


NOo 


NO  + 0 - 2.60  eV. 


as  a function  of 


The  momentum- transfer  cross  section  for  a (12-6-4)  potential  has 

2 

been  tabulated  in  the  dimensionless  form  q,,  (v)/irr 
12  m 

j p\'  /e  and  the  parameter  y (References  17,  18).  Figures  4-2  to  4-8 

show  plots  of  the  cross  sections  that  are  found  by  interpolation  from 

In  Section  5 the  rate  coefficients  s^(Tm>V)  and 


h.,  (T  ,V) 
jk  v m’  1 


this  table. 

for  elastic  scattering  are  found  by  integration  from  Eqs.  2-33  and  2-34; 

they  are  plotted  in  Figures  5-1  and  5-3.  An  alternate  method,  which  we  used 

as  a check, is  to  find  s (T  ,0)  by  interpolation  from  the  tabulated  collision 
— fill  m 

integral  ’ 1 (T)  in  Reference  17  together  with  the  relation  (see  Refer- 
ence 2,  p.139) 


sjk(T,0) 


JL 

3/i FUVV/ 


i) 


(T) 


(3-21) 


For  0,  + 0,  for  which  the  nature  of  the  interaction  is  less 

certain,  the  estimated  uncertainty  in  s (or  K^)  is  a factor  of  two 
(Reference  14)  for  temperatures  up  to  0.22  eV  and  an  order  of  magnitude 
at  high  temperatures.  For  N+  + and  Of  * 0^  we  feel  that  the  un- 
certainty should  be  a factor  of  two  over  the  whole  temperature  range,  and 
for  the  first  five  systems  in  Table  3-4  we  estimate  ±30%  uncertainty. 


4! 


SS^SSSHS©}^ ^^mS^S^i^^s^s^iSKm'isss 


(c)  Scattering  from  50  to  30000  eV 

In  ion-atom  encounters  at  ccnter-of-mass  energies  above  200  eV 
elastic  scattering  is  exceeded  by  charge  exchange  (see  Figures  4-1  to  4-8). 
These  charge-exchange  encounters  cause  very  little  momentum  transfer  between 
the  fast  particles  and  the  slow  ones.  The  slowing  down  of  the  fast  particles 
is  most  likely  brought,  about  by  atom-atom  scattering.  We  have  therefore  con- 
centrated on  atom-atom  scattering  rather  than  ion-atom  scattering  in  this 
energy  range. 

We  use  a two-atom  Thomas -Fermi  interaction  potential  with  a scaling 
rule  for  unlike  atoms  (Reference  19).  The  range  of  validity  of  the  treatment 
is  restricted  at  the  lower  end  to  energies  much  gi'catcr  than  molecular  binding 
energies,  say  center-of-mass  kinetic  energies  above  50  eV.  At  the  upper  end 
inelastic  scattering  involving  excitation  and  ionization  overtakes  elastic 

g 

scattering  at  a speed  around  10  cm/sec  (see  Appendix  C and  Figure  4-6). 


The  colliding  pair  is  described  simply  by  the  atomic  numbers  Z^ 
and  Z9.  The  scales  for  distance  and  energy  are 


1_  /3ir\  2 ' '~7  2 \2  / 3 

wv^)  a° 

0.469  x 10'8  / — \ 

\>lh  + vv 


(3-22) 


eo  - e'Z^/A 


(3-23) 


= 0.0307  Z, 
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The  scaled  centcr-of-mass  kinetic  energy  of  the  colliding  particles  is 
denoted  by 

W = ~2  I'V2/^  , (3-24) 

and  the  momentum- transfer  cross  section  is  given  in  the  dimensionless  form 
l;(W)  in  Table  .3-5,  which  is  taken  from  Table  2 of  Reference  19: 

q„(v)  = X21(W).  (3-25) 

At  intermediate  energies  we  use  log-log  Interpolation  in  the  table,  and  at 
the  high-energy  end  we  extrapolate  by 

F(W)  = 0.066(10/W)5/2  if  W > 10.  (3-26) 

These  cross  sections  for  II  I + II  , N + N , 0 + 0,  and  N + 0 are 

c c 

shown  in  Figures  4-1,  4-2,  4-3,  and  4-5,  respectively;  these  are  believed 
to  be  accurate  to  better  than  a factor  2 (Reference  19).  For  atom-molecule 
collisions  such  as  0 + 0^  we  have  simply  added  the  appropriate  atom-atom 
scattering  cross  sections  (Figures  4-6  to  4-8). 


Collisions  of  singly  charged  ions  with  atoms  are  more  complicated 
than  atom-atom  collisions  because  one  should  take  into  account  charge 
exchange.  This  is  not  done  in  the  classical  scattering  calculation  used 
here  tor  atom- atom  cross  sections.  For  a rough  estimate  one  might  take 
the  ion-neutral  scattering  cross  section  to  equal  the  corresponding  neutral - 
neutral  one;  we  do  this  in  Section  5 in  extending  the  rate  coefficients 
and  hjk  to  temperatures  above  50  eV. 

Table  3-5.  The  function  F(W)  = qD(v)/A2  (Reference  19,  Table  2). 


w 

Fj-WJ .... 

1 

W 

nmn 

6.005 

no. 

i 

0.3 

5.28 

0.010 

70. 

1. 

1 .46 

0.025 

37. 

3. 

0.367 

0.05 

22.6 

10. 

0.066 

JL1 

J3.4_ 



>in 

£LH66*|liUHJ^2 
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1.1 


SECTION  4 

CHARGE-EXCHANGE  CROSS  SECTIONS  FOR  ATMOSPHERIC  SPECIES 


INTRODUCTION 


We  consider  symmetric  charge  exchange  of  six  positive  ions: 


He  + lie  ■*■  He  + He  , (4-1) 


N ■*■  N ->-  N + N , (4-2) 


0+  + 0 -►  0 + 0+  , (4-3) 


N2+  + N2  N2  + N2+  » (4-4) 

N0+  + NO  -*•  NO  + N0+  , (4-5) 


°2+  + °2  ■*■  °2  + °2+  * (4-6) 


4 

a 

_! 

.a 

$ 


;wi 


M. 

m 

• — -g 

*v 


x-y 

M 


$ 


■m 


1 


and  eight  asymmetric  charge-exchange  reactions: 


N+  + 0 -*■  N - 0+  + 0.92  eV  , 


+ N2  + N + N0  - 1.05  eV  , 


+ N ■*  N2  + N + 1.05  eV  , 


3+  (4s)  + 0,  •>  0 + 0,+  + 1.56  eV  , 


(4-7) 

(4-8) 

(4-9) 

(4-10) 
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o'  i2n)  + 02  -►  0 + 02+  4-  4.88  cV  , 


(4-11) 


0*  + 0 -►  02  + o'  - 1 .56  cV  , 


(4-12) 


i\''  + NO  -►  N + N0+  + 5.27  eV  , 


O'  (4S  ) + NO  -►  0 + N0+  + 4.35  oV  , 


(4-13) 


(4-14) 


o'(2ll)  + NO  > 0 + N0+  + 7.68  eV  , 


(4-15) 


N+  + U ->  N + U+  + 8.4  eV  , 


O'  + U * 0 + U*  • 7.5  CV  . 


(4-16) 


(4-17) 


For  three  other  asymmetric  reactions  we  have  searched  the  literature  but 
found  no  reliable  cross  section: 


o'  4-  N -►  0 4-  N+  - 0.92  eV  , 


(4-18) 


NO'  + N NO  + N+  - 5.27  eV  , 


(4-19) 


NO  4-  0 -*■  NO  + 0 - 4.35  eV  . 


(4-20) 


Ground-state  reactants  are  emphasized;  information  on  the  states  of  the  re- 
actants is  included  where  available.  The  products  of  the  asymmetric  re- 
actions are  likely  to  be  excited,  and  all  states  are  included.  The  energy 
excesses  shown  are  only  limits  because  they  are  calculated  by  assuming  that 
the  products  as  well  as  the  reactants  are  in  their  ground  states. 


We  have  searched  the  literature  through  November  1974  for  measured 
cross  sections.  Where  exp  . imenta.l  data  is  available  we  have  relied  on  it 
rather  than  on  theoretical  calculations. 
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Tables  4- la  to  4- lk  list  the  studies  selected.  The  first  column 
includes  the  name  of  the  first  author  and  the  publication  date.  The  second 
column  gives  the  experimental  technique  and  the  particle  detection  method, 
whore  important,  as  follows: 

IA  flowing  afterglow 
SA  stationary  afterglow 
l)T  drift  tube 

I'DT  flowing-afterglow  drift  tube 
ICU  ion  cyclotron  resonance 

UPMS  high-pressure  mass  spectrometer  ion  source,  operated 
in  the  many-col lision  pressure  regime 
15  beam  in  gas  cell,  operated  in  the  single-collision 
pressure  regime 

liq  beam  in  gas  cell,  equilibrium  method,  fast-particle 
detection 
Cl!  crossed  beams 

MB  merging  beams,  including  magnetic,  analysis  of  products 
f detection  of  fast  particles 

s detection  of  slow  particles.  If  the  energy  of  the 

5 

primary  ions  in  a beam  experiment  exceeds  10'  eV, 
t ho  slow  product  ions  include  many  that  are  produced 
by  ion-impact  ionization.  Also,  if  the  energy  of  the 
primary  ions  is  very  low,  the  slow  ions  include 
elastically  scattered  primary  ions, 
m product  ions  analyzed  by  a mass  spectrometer 

i reactants  isotopically  labeled. 


Another  column  gives  the  range  of  ion  energy  for  which  a cross  section  is 

2 

measured.  If  the  author  used  ccntcr-of-mass  energy  E^  = (m^m^/On  ■*-m.,))vi72, 
this  range  is  given  in  brackets.  Where  a temperature  range  is  given,  the 
measured  quantity  is  indicated  under  Remarks;  it  is  cither  the  ion  mobility 


(see  liq.  3-20)  in  the  case  of  ions  in  their  parent  gas  or  the  rate  coef- 
ficient in  the  case  of  asymmetric  charge  exchange. 


4.2 


SYMMETRIC  CHARGE  EXCHANGE 


■V  4 

He  + He  •>  He  * lie  : We  selected  the  investigations  in  Table 

4-  la.  Reference  26  is  a study  of  3lle+  + 'lie  •►'Hie  + Hlc^,  which  is  not 
strictly  symmetric,  but  the  difference  between  the  ionization  energies  of 
3llc  and  'lie  is  assumed  to  he  too  small  to  affect  the  resonance  character  of 
the  charge  exchange.  References  28  to  40  agree  well  and  provide  a cross 
section  down  to  ion  energy  4 eV,  which  is  shown  in  figure  4-1.  For  con- 
venience Figures  4-1  to  4-8  are  placed  together  at  the  end  of  Section  4, 
immediately  preceded  by  the  notes  on  these  figures. 


3 3 

lor  1.1  x 10'  < v < 2.3  x 10'  cm/sec  we  use  the  mobility 
measurements  of  Reference  25  together  with  F.qs.  3-20,  2-33b,  and  2-43  to 
obtain  an  empirical  power-law  fit  to*  clp(v)  + 2 o^lv): 


q0(v)  + 2 o ( v ) - 92  x io‘16  cm2  (v/1.58  x io5  cm/sec)'0'56  . 

(4-21) 

At  those  low  energies  there  is  no  experimental  data  on  how  to  separate 
tl\is  into  q (v)  and  o ^ ( v ) individually;  we  separate  them  arbitrarily 
by  taking  the  elastic  scattering  q^(v)  to  be  given  by  its  polarization 
limit,  liq.  3-15.  For  ion  energies  below  1 cV  this  charge-exchange  cross 
section  °x(v)  is  dotted  in  Figure  4-1  as  a reminder  that  it  was  chosen 
merely  to  make  the  total,  q^(v)  + 2 o^(v),  consistent  with  measurements. 


Strictly  one  obtains  a fit  to  the  quantity 

do 


/uo 

(l  + cost))  — dS) 


4 it 


,(1 


which  is  the  Ql  } of  Reference  J,  S-7-l, 
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Table  4-la. 


Studies  of  the  cross  section  for  He+  + He  He  + He+. 


Reference 

Technique 

Ion  Energies 

Remarks  1 

23  (Heiche,  1970) 

theory 

l2xl°'5<Ean<4  eV) 

24  (Dickinson,  1968) 

theory 

l0.0037<Ecn<0.13  eV] 

25  (Orient,  1967) 

DT 

1190°<T<650°K) 

Ion  Mobility 

26  (Smith*  1974) 

ICR.i 

[0.1<Ecm<20  eV) 

Used3lle+  +4lle 

27  (Mahadevan,  1968) 

B,s 

1 to  75  eV 

28  (Cramer,  1957) 

B,s 

4 to  400  eV 

29  (Belyaev,  1968) 

MB 

L7<Ec,n<10°  eV] 

30  (Hayden,  1964) 

B,s 

50  to  1000  eV 

31  (Ghosh,  1957) 

B > S 

100  to  8C0  eV 

32  (Nagy,  1969) 

B,f 

400  to  2000  eV 

33  (Shelton,  1971) 

B,S 

2 to  20  keV 

34  (Stedeford,  1955) 

B,s 

2 to  40  keV 

35  (Barnett,  1958) 

Eq 

8 to  200  keV 

36  (DeHeer,  1966) 

B,s 

10  to  140  keV 

37  (Fedorenko,  1956) 

B,s 

13  to  180  keV 

38  (Jones,  1959) 

B,f  ,nt 

25  to  100  keV 

39  (Allison,  1956) 

B,f 

100  to  450  keV 

40  (Nikolaev,  1961) 

B,f  ,m 

320  to  1300  keV 

a 


I 

•s 


I 

i- 

4 

i 

I 

-T 

ji 

"1? 


V3i 
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+_N  N + N* : We  use  the  only  known  measurement  (Reference  25)) 

and  extrapolate  in  both  directions  by  following  the  shape  of  the  Rapp  and 
Francis  theoretical  curve  (Reference  II) . See  Table  4-ll>  and  Figure  4-2. 
Above  30  kcV  ion  energy  wc  estimate  the  cross  section  from  measurements  by 
Fite's  group  on  the  charge  exchange  of  N + N0  and  0*  « 0,  (Reference  73) 
and  of  0+  * 0 (Reference  44);  at  a given  speed  our  estimate  is 

4 + V°++0) 

o(N  +-N)  * o (N  vN?)  — — — . (4-22) 

“ 0(0  +0,) 


Table  4-lb.  Studies  of  the  cross  section  for  + N N + Nf 


Reference 

Technique 

Energy  Range 

29  (Belyaev,  196H) 
41  (Rapp,  1962) 

MB 

Theory 

7 < EC||1  < 100  eV 
0.005  < E < 35000  eV 

0*  4-0  -»  o r o* ; The  agreement  among  the  measurements  in  Table 
lie  is  very  good,  and  the  result  is  shown  in  Figure  4-3.  We  have  again 
extrapolated  down  to  0.3  eV  ion  energy  by  following  the  shape  of  the 
theoretical  cross  section  of  Reference  II. 

Oxygen  atoms  formed  by  thermal  dissociation  of  0o  are  in  the 
ground  state,  while  those  produced  by  an  rf  discharge  are  likely  to  include 
excited  atoms.  For  the  resonant  reaction  1-3  one  would  not  expect  the 
presence  of  excited  states  to  influence  the  cross  section  greatly 
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Table  4-lc.  Studies  of  the  cross  section  for  0*  + 0 -1-  0 + 0 . 


Re ference 

Technique 

Production 
of  0 atoms 

Ion  Energies 

42  (Stebbincjs,  1%4) 

611,5,111 

rf  discharge 

40  to  10000  eV 

43  (Rutherford,  1974) 

C(l,s,m 

thermal  dissoc. 

60  to  500  eV 

44  (Lo,  1969) 

B,t,m 

thermal  dissoc. 

30  to  2000  keV 

N.,  + N,  N?  + N,  : We  selected  fifteen  sources  of  experi- 

mental cross  sections,  as  listed  in  Table  4-ld.  Discrepancies  among 
those  cross  sections  are  partly  due  to  differences  in  the  distribution 
of  electronic  and  vibrational  states  of  the  primary  N,+.  In  all  of  the 
beam  experiments  in  this  table  except  Reference  31  the  primary  N.,+  ions 
were  produced  by  electron  impact  on  the  electron  energy  is  in  the 

third  column,  lour  investigators  (References  52,  54,  57  and  58)  studied 
the  effect  of  changing  the  distribution  of  states  of  the  primary  N.,+  by 
varying  this  electron  energy.  They  find  that  changing  the  electron  energy 
from  1G  or  18  c V to  high  energies  (25  to  DO  eV)  decreases  the  cross  section 
somewhat.  The  largest,  decrease  measured  was  a factor  1.35. 


Maier  (Reference  52)  has  calculated  the  initial  distribution  of 

electronic  and  vibrational  states  of  N,  for  three  electron- impact 

+ y 

energies.  Ot  the  known  electronic  states  only  N,  (A"ll  ) is  sufficiently 


populated  and  has  a long  enough  lifetime  to  affect  the  measured  cross 

+ y 

sections.  The  cross  sections  tor  N.,  (A  'll  ) ions  in  various  vibrational 

levels  to  charge-transfer  with  ground-state  N , are  estimated  theoretically 

t.o  bo  significantly  smaller  than  the  cross  sections  for  ground-state 
+ > ( 

N,  (“■•„)  reacting  with  ground-state  N.,  (Reference  58,  p.  SSOG).  The 
- »*>  •-  > 
lifetimes  of  the  vibrational  levels  of  the  A"n  state  (Reference  52, 

Table  11)  are  comparable  to  the  elapsed  time  between  formation  and  reaction 


Table  4-ld.  Studies  of  the  cross  section  for  N./  + "•  N,,  + N^+. 


Reference 

Technique 

Electron 

Energy 

Ion  Energies 

Remarks 

45  (Moseley,  1969) 

il  l ,111 

[T=300"Kj 

Ion  Mobil i ty 

46  (Huntress,  1971) 

ICR 

[300‘VT.<50QoK] 

Note  b 

47  (Kobayashi,  1975) 

iii.in 

f0.05<E  <2.5  eV] 

46  (Ticrnan) 

B,i 

not  known 

-0.3  eV 

49  (Nichols,  1966) 

B,s 

96  eV 

0.5  to  17  eV 

50  (Lovonthal , 1967) 

U,s,m,i 

60  to  70  eV 

1 to  50  eV 

61  (St.ehhings,  1963) 

CB.s 

•90  eV 

15  to  5000  eV 

52  (Maicr,  19/4) 

D.s.m.i 

19,2  eVa 

20  to  45  oV 

53  (Gustafsson,  1960) 

B,s,m 

high,  -70  eV 

50  to  900  eV 

54  (Amine,  1964) 

13,  s 

21 . 5 eVa 

50  to  1000  0 V 

31  (Ghosh,  1957) 

l),s 

rf  discharge 

100  to  «00  eV 

55  (Savage,  I960) 

MB 

30  to  50  eV 

100  to  1000  eV 

56  (Neff,  1964) 

13, s 

not  stated 

190  to  3000  eV 

57  (McGowan,  1964) 

13,  f 

55  eVa 

600  to  2000  0 V 

50  (Flannery,  1973) 

- 

B.f 

Note  a 

760  to  2200  eV 

Note  c 

ll  The  electron  energy  was  varied  to  study  the  effect  of  changes  in  the  distribution 
of  states  of  the  primary  ions, 

The  quantity  plotted  in  Huntress's  figure  3 is  k = :'/n,  where  is  the  velocity 
average  of  the  (momentum-transfer)  collision  frequency  and  n is  the  number 
density  of  N?.  This  is  related  to  the  reduced  mobility  K0  by  Huntress’s 
Eg.  24,  and  our  Eg.  3-20  then  gives  Sj^  + si;*  = m^  k.  The  abscissa  of  Mg.  3 
is  the  ion  kinetic  temperature  Tj;  the  temperature  of  the  neutrals  wo  take  to 
be  T„  293  °K. 

Gives  only  relative  cross  sections  as  a function  of  electron  energy. 

T1 


■51 


is 


c 


of  the  N + in  various  beam  experiments.  This  latter  time  depends  on  ion 
energy  and  is  not  stated  for  most  of  the  beam  experiments  in  Table  4- Id. 

Thus  one  can  find  only  an  upper  limit  to  the  fractional  population  of  the 

2 . + 

A n state  in  the  incident  N.  beam, 
u 2 

We  have  not  attempted  to  adjust  the  measured  cross  sections  to 
refer  to  a given  distribution  of  states,  such  as  a pure  beam  of  ground- 
state  N_+  ions.  The  discrepancies  due  to  different  distributions  of 
states  in  the  primary  N seem  to  be  no  larger  than  the  discrepancies 
due  to  unknown  other  sources.  For  example,  at  an  ion  energy  around  50  eV 
the  cross  sections  ft'om  Table  4-ld  differ  by  factors  up  to  2.0,  and  the 
smallest  (Reference  52)  and  largest  (Reference  54)  are  from  experiments 
that  used  similar  low  electron  energy. 

Figure  4-4  shows  our  recommended  cross  section.  The  extension 
to  ion  energies  above  5000  eV  is  an  extrapolation. 

4 7 

For  5.9  x 10  < v < 7 < 10  cm/scc  we  use  the  dc  and  ac  mobility 
measurements  of  References  45  and  46  together  with  Fqs.  5-20,  2-33b,  2-21, 
and  2-43  to  obtain  an  empirical  power- law  fit  to  q^(v)  + - (v) : 

q()(v)  + 2 Ox(v)  = 1.6  x io  cm^  (v/5.97  * 10^  cm/sec)  ^ ' * 

(4-23) 

We  calculated  0 below  1 eV  by  the  same  arbitrary  separation  procedure 
x + + 

that  is  described  above  for  lie  + lie  He  + He  , and  the  same  remarks 
apply. 


0 + + 0^  ► 0 + 0 +:  The  seven  references  in  Table  4-le  show 

- c t + 

>■  ? serious  discrepancies.  The  07  + 0?  measurements  of  Reference  47 

were  published  after  Figure  4-4  was  completed,  but  only  small  adjustments 
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Table  4-le,  Studies  of  the  cross  section  for  Og  + 0^  -*■  0^  + . 


Reference 

Technique 

Ion  Energies 

Remarks 

59  (Snuggs,  1971) 

0T  ,m 

[T  = 300  °K] 

Ion  Mobility 

47  (Kobayashi,  1975) 

DT  ,m 

[° . 04 5 < Ecm  < 2.8  eV] 

60  (Varney,  1970) 

DT  ,m 

-0.2  eV 

48  (Tiernan) 

B,i 

-0.3  eV 

51  (Stebbinas,  1963) 

CB,s 

25  to  6500  eV 

61  (Ghosh,  1957) 

B,s 

100  to  800  eV 

54  (Amine,  1964) 

B , s 

100  to  1000  eV 

5 

are  indicated  for  speeds  above  2 * 10  cra/sec.  Below  1 eV  ion  energy  we 
used  the  mobility  measurement  of  Reference  59  and  followed  the  same 
procedure  described  for  lie*  + He  amd  + N,,.  The  result  is  more  un- 

certain in  this  case  because  the  temperature  dependence  of  the  mobility  was 
not  measured  and  must  be  estimated. 

Two  investigators  (References  62  and  54)  have  studied  the  effect 
of  changing  the  distribution  of  states  of  the  primary  0,+  by  varying  the 
electron  energy  used  in  producing  the  0 + by  electron-impact  ionization. 
Compared  to  the  accuracy  of  the  measurements,  this  effect  is  unimportant. 

NO*  + NO  ->  NO  + NO* : The  cross  sections  from  Table  4- If  are 

considerably  more  uncertain  than  for  the  other  cases  of  symmetric  charge 
exchange  considered  in  this  study.  One  cause  is  the  presence  of  excited 
states  in  the  primary  NO+  ions  (References  51  and  62).  For  primary  NO 
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Table  4-1  f.  Studies  of  the  cross  section  for  N0+  + NO  ■*■  NO  + N0+. 


Reference 

Technique 

Ion  Energies 

Remarks 

63  (Volz,  1971) 

DT,m 

[T  = 300  °K] 

Ion  Mobility 

64  (Henglein,  1962) 

B,s  ,111 

5-40  eV 

51  (Stebbings,  1963) 

CB , s ,111 

25  to  4000  eV 

61  (Ghosh,  1957) 

B,s 

100  to  800  eV 

62  (Moran,  1974) 

B,f 

400  to  1200  eV 

Note  a 

a Gives  only  relative  cross  sections  as  a function  of  electron  energy. 


ions  produced  by  electron-impact  ionization,  varying  the  electron  energy 
from  9 to  above  30  eV  increases  the  cross  section  by  a factor  1.5  at  ion 
energy  400  eV  (Reference  62);  this  effect  is  opposite  to  that  for  N?  + N? 
mentioned  above. 

4.3  ASYMMETRIC  CHARGE  EXCHANGE 

N*  + 0 -»  N + 0* : from  the  sources  in  Table  4-lg  we  put  together 

the  cross  section  in  Figure  4-5.  The  two  arrows  above  this  curve  are 
measured  upper  limits.  In  References  43  and  67  the  0 ->101113  are  known  to 

be  in  the  ground  state  because  they  are  formed  by  thermal  dissociation.  The 
cross  section  for  reaction  with  the  metastable  N+(*D)  ion, 

nVd)  + 0 -►  N + 0+  , (4-24) 

+ 

is  reported  to  be  not  large  compared  to  that  for  ground-state  N ions, 
even  at  ion  energies  as  low  as  5 eV  (Reference  43) . 
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Table  4- 1 g . Studies  of  the  cross  section  for  N + 0 -*•  N + 0 . 


Reference 

Technique 

Production  of 
0 atoms 

Ion  Energies 

Remarks 

65  (Bortner,  1973) 

not  stated 

not  stated  d 

Estimate  of 
rate  coeff. 

66  (Neynaber) 

MB 

charge  transfer 

0.19  eV 

Upper  1 imit, 

ox<10"17cm2 

( 5 eV 

Upper  limit 

43  (Rutherford,  1974) 

CB,s,m 

thermal  dissoc. 

(400  eV 

67  (Lo,  1970) 

B,f ,m? 

thermal  dissoc. 

30  to  2000  keV 

a We  have  assumed  that  the  estimate  applies  for  200  ° < T < 700  °K. 


N*  + U ■»  N + U+  and  0*  * 0 ->  O + IJ* : Also  shown  in  Figure  4-5 

are  1974  results  of  crossed-beam  experiments  of  Neynaber,  Rutherford,  and 
Vroom  (Reference  68). 

N+  + -*■  N + N7+:  1,1  the  detailed  low-energy  study  of  Maier  and 

Murad  (Reference  69)  the  reactants  were  isotopically  labeled  and  the  kinetic 
energies  of  the  product  N^,+  were  determined.  They  conclude  that  this 
reaction  proceeds  by  both  atom  and  electron  transfer.  For  ion  energies 
below  16  eV  the  cross  sections  for  electron  and  atom  transfer  are  roughly 
equal,  but  at  higher  energy  the  cross  section  for  electron  transfer  is 
dominant.  Also  the  kinetic  energy  transferred  into  internal  energy  of  the 
products  is  considerably  larger  than  the  minimum  of  1.05  eV  required  for 
this  reaction  to  occur;  it  increases  from  3.2  to  10  eV  as  the  ccnter-of- 
raass  kinetic  energy  increases  from  4 to  50  eV. 
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The  effect  of  excited  states  of  N+  has  been  studied  in  References 
71  and  53.  Neynaber  et  al.  (Reference  71)  conclude  that  the  cross  section 
for  the  reaction 

N+(1D)  + N2  ->  N * N2+  (4-25) 

is  expected  to  be,  and  is  observed  to  be,  small.  (They  used  ion  energies 
of  ~10  to  500  eV.)  McGowan  and  Kerwin  (Reference  53)  report  only  a small 
increase  (-20  percent)  in  the  cross  section  as  the  electron  energy  used  to 
produce  the  N ions  is  varied. 

The  differences  among  the  three  cross  sections  presented  in 
References  69  and  70  are  a factor  of  2 to  4,  and  when  those  of  References  51 
and  53  are  considered,  the  discrepancies  are  far  larger.  See  Table  4-lh. 
Maier  and  Murad  (Reference  69)  discuss  the  differences  and  recommend  that 
the  cross  section  of  Reference  70  be  preferred  when  the  primary  N+  ions 
are  produced  by  electron  impact  on  low-pressure  N^.  We  adopt  the  cross 
section  of  Reference  70  for  ion  energies  up  to  500  eV. 

At  high  energies  the  agreement  between  References  72,  73  and  40 

is  good.  (Note  that  in  References  40  and  72  the  cross  sections  are  given 
2 2 

in  cm  /atom  instead  of  cm  /molecule  as  used  in  Figure  4-6).  For  these 
energies  the  cross  sections  refer  to  the  reaction 

N+  + N?  N + slow  products  . (4-26) 

Figure  4-6  shows  the  recommended  cross  section,  in  which  these  high-energy 
results  arc  joined  to  that  of  Reference  70  by  a smooth  dotted  curve. 
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Table  4-1 h.  Studies  of  the  cross  section  for  N+  + N2  + N + N?  . 


Reference 

Technique 

69  (Maier,  1971) 

( B,s,m,i 
| B , s ,m 

70  (Neynaber,  1971) 

CB,s,m 

51  (Stebbings,  1963) 

CB,S 

53  (Gustafsson,  1960) 

B,s,m 

57  (McGowan,  1964) 

B,f 

72  (Ormrod,  1971) 

B,f  ,m 

73  (Brackmann,  1968) 

B,f  ,IH 

40  (Nikolaev,  1961 ) 

B , f ,m 

Ion  Energies 

3 to  50  eV 
1.5  to  100  eV 
1.5  to  500  eV 
20  to  5000  eV 
25  to  900  eV 
2000  eV 
15  to  70  keV 
30  to  1300  keV 
500  to  2400  keV 


N + + N -*-  ♦ N+ : Kretschmer  and  Petersen  (Reference  74)  have 

measured  the  rate  coefficient  at  300  ’K  temperature  by  an  indirect  technique 

- 13  3 

using  Langmuir  probes.  Their  rate  coefficient,  1.4  x 10  cm  /sec,  is 
much  lower  than  the  upper  limit  determined  by  the  afterglow  experiments  of 
Ferguson  et  al.  (Reference  75).  The  result  of  Reference  74  is  shown  as  a 
point  on  Figure  4-6. 

0+  + 0>  -►  O + 02+:  The  cro,->c  sections  for  ground-state  O (4S) 

and  for  met  as  t. Me  0+  (V)  have  been  measured  separately  in  three  beam 
experiments  (References  48,  68,  and  85;  see  Table  1-li).  Their  results  are 
consistent,  and  we  have  shown  a separate  cross  section  for  0 (2D)  + 02 
in  Figure  4-7. 
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Table  4-li.  Studies  of  the  cross 

section  for  0+  + 02 

- o + o2+. 

- A 

3 

Reference 

Technique 

Ion  Energies 

Remarks 

'?£ 

76  (Dunkin,  1968)  ) 
| ( 

[80°<T<600°K] 

FA 

Rate  Coefficient 

:3f, 

77  (Ferguson,  1969)) 

M 

78  (Smith,  1968)  ) 

1 

( 

SA 

[180°<T<600°K] 

Rate  Coefficient 

i 

79  (Copsey,  1966)) 

I 

80  (Lindinger,  1974) 

FA 

[300°<T<800°K] 

Rate  Coefficient 

..-6 

81  (Warneck,  1967) 

HPMS 

[700°<T.<1100°K]a 

Rate  Coefficient 

J 

82  (Pharo,  1971) 

ionospheric 

[T  = 900° K] 

Rate  Coefficient 

83  (Johnsen,  1973) 

DT  ,m 

300° K to  2 eV 

:yd 

-I* 

•S' V 

84  (McFarland,  1973) 

FDT  ,m 

[0.05<Ecm<3.2  eV] 

48  (Tiernan) 

B 

~ 0.3  eV 

Also  0+(2D)+02 

■ 

68  (Neynaber,  1974) 

CB,s  ,m 

4 to  500  eV 

Also  0+(2D)+02 

M-C* 

A 

85  (Stebbings,  1966) 

CB,s,m 

10  to  100  eV 

Also  0+*+02 

: S 

57  (McGowan,  1964) 

B»f 

1.6  keV 

: A^i 

; M 

86  (Turner,  1969) 

B 

8 to  60  keV 

87  (Solov'ev,  1972) 

B,f,m 

28  to  180  keV 

: f 

73  (Brackmann,  1968) 

B,f  ,m 

30  to  2000  keV 

■t 

a To  compare  to  other  temperatures  in  this  column  one  must  calculate 
the  temperature  Tm  (Eq.  2-21  or  2-47)  characterizing  the  distribution 
of  relative  velocities  of  0+  - 0?  pairs.  For  an  0?  temperature  of 
300°K  the  range  of  this  temperature  is  600°<Tm<80(r K. 
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The  drift- tube  measurements  are  believed  to  refer  to  ground-state 
0+  (*S)  ions  (Reference  83) , and  we  have  assumed  the  same  for  the  afterglow 
experiments.  The  temperature-dependent  rate  coefficients  have  been  used 
to  estimate  charge- exchange  cross  sections  °x(v)  by  using  Eq.  B-llb  to 
obtain  an  empirical  power- law  fit  to  cr^fv).  Some  of  the  recent  measure- 
ments (References  80,  82,  and  48)  indicate  a lower  ci'oss  section  at  low 

energies  than  do  References  76  to  79,  81,  and  83;  the  discrepancy  is  nearly 

5 8 

a factor  of  2 in  the  range  1 * 10  < v < 2 * 10'  cm/sec.  We  choose  an 
intermediate  cross  section  here  and  indicate  the  conflict  by  a dashed  line 
in  Figure  4-7, 

At  intermediate  energies  we  follow  the  cross  section  of  Noynaber 
et  al.  (Reference  68).  At  high  energies  References  87  and  73  agree  well 
and  we  have  followed  the  latter.  For  these  energies  the  cross  sections 
refer  to  the  reaction 

0*  + 0 + slow  products  . (4-27) 

In  Figure  4-7  this  high-energy  curve  is  connected  to  that  of  Reference  68 
by  a smooth  dotted  curve. 

0?+  + 0 0+  + 02:  Hie  only  measurement  is  the  crossed-beams 

experiment  of  Stebbings  ct  al.  (Reference  88),  covering  ion  energies  from 
70  to  6000 eV.  This  measurement  is  probably  strongly  affected  by  the  presence 
of  excited  atoms  and  excited  ions  in  the  incident  beams.  To  estimate  the 
cross  section  for  ground-state  reactants  we  compared  measurements  of 
Reference  88  to  those  of  Reference  43  for  two  similar  charge-exchange 
reactions,  those  for  N^+  + 0 and  N+  + 0 (reaction  4-7).  In  Reference  43, 
where  the  reactant  O atoms  are  in  the  ground  state  and  the  ions  are  pro- 
duced by  impact  of  lower  energy  electrons,  these  two  cross  sections  are 
smaller  than  those  of  Reference  88  by  a factor  of  about  6.  For  our  estimate 
of  the  cross  section  for  0^+  + O,  wc  divided  the  cross  section  of  Reference 
88  by  6 (and  at  2 keV  ion  energy  by  a more  conservative  factor  of  3).  Figure 

7 

4-7  shows  the  result,  which  covers  only  the  range  v < 2 * 10  cm/sec. 
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N*  -*•  NO  -*■  N ■*•  NO+:  Rutherford  (Reference  91)  has  extended  the 


crossed- beans  measurements  of  Turner  (Reference  92)  down  to  1 cV  ion  energy; 

At  room 


“1C  2 

the  ci*oss  section  remains  nearly  constant  at  a **  1.5  * 10*  1 enr 


temperature,  however,  the  flowing  afterglow  measurement  (Reference  89)  in- 
dicates a cross  section  of  about  100  * 10 cm‘  for  v *s  7 x 10^  cm/sec. 
Between  room  temperature  and  1 eV  are  the  measurements  of  References  90  and 
48  at  two  well-placed  energies.  It  is  unfortunate  that  the  range  of  ion 
temperature  or  mean  ion  energy  used  for  this  reaction  by  Wameck  (Reference 
90)  is  not  stated  more  precisely,  but  his  measurement  is  nonetheless  signi- 
ficant. The  cross  sections  from  the  four  experiments  mentioned  so  far  (the 
first  four  in  Table  4-1 j)  decrease  monotonicaliy  as  the  energy  is  increased 
to  1 cV.  Taken  together,  they  support  the  steep  drop  in  cross  section 
shown  in  Figure  4-8. 


It  is  not  known  whether  atom  transfer,  as  opposed  to  election 
transfer,  is  important  at  low  energies  (Reference  89).  The  effect  of 
excited  states  of  N+  was  studied  in  Reference  92  by  varying  the  energy 
of  the  electrons  used  to  produce  the  N+;  the  cross  section  is  relatively 
insensitive  to  electron  energy. 


At  high  energies  we  average  the  charge-exchange  cross  sections 

of  Reference  73  for  N+  + N_  (reaction  4-8)  and  for  N+  + 0-.  Figure  4-8 

L 1 4 8 

shows  our  recommended  cross  section,  which  covers  the  range  t>  * 10  < v < 5 * 10 

cm/sec.  (At.  high  energies  this  curve  is  indistinguishable  from  that  for 

0+  +■  NO). 


0+  + NO  -*•  0 + N0+ : 


Both  beam  expci iments  in  Table  4- lk  have 

+ ,4C 


measured  separately  the  cross  sections  for  ground-state  0 ( S)  and  for 
metastable  0+(2L)),  and  the  results  of  Reference  48  are  consistent  with 
extrapolations  of  the  results  of  Reference  68.  At  the  lowest  ion  energies 
in  these  experiments  (0.3  and  1.5  eV) , however,  even  the  measurements  for 
ground-state  0+  exceed  the  flow-drift  tube  result  (Reference  93)  by  a 
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Table  4-1 j.  Studies  of  the  cross  section  for  N+  + NO  -►  N + N0+. 


Reference 

Technique 

Ion  Energies 

Remarks 

89  (Goldan , 1966) 

FA 

[T«300°K] 

Rate  Coefficient 

90  (Warneck,  1967) 

HPMS  or  LPMS 

Note  a 

Rate  Coefficient 

48  (Tiernan) 

B,i 

~0.3  eV 

91  (Rutherford) 

CB,s  ,m 

1 to  > 5 eV 

92  (Turner,  1966) 

CB,s  ,m 

5 to  200  eV 

a 

The  ion  temperature  range  for  the  group  of  rate  coefficients  sum- 
marized in  Ref.  90  is  400°<T|<2000°K,  but  each  individual  reaction 
was  studied  only  over  part  of  this  range.  See  footnote  to  Table  4-li. 


Table  4-1  k.  Studies  of  the  cross  section  for  0+  + NO  -*•  0 + N0+. 


Reference 

Technique 

Ion  Energies 

Remarks 

93  (McFarland,  1974) 

FDT  ,ri 

| [T=280°K] 

Upper  limit  for 
rate  coefficient 

48  (Tiernan) 

B 

|[0.22<Ecm<4  eV] 
-0.3  eV 

Also  0+(2D)+N0 

68  (Neynaber,  1974) 

CB,s,m 

1.5  to  500  eV 

Also  0+(2D)+N0 
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factor  of  3 or  more.  Being  unable  to  resolve  this  discrepancy  we  give  an 
intermediate  cross  section  for  0*('*S)  + 0^  at  these  energies  and  indicate 
the  decreased  accuracy  by  a dotted  line  in  figure  4-8.  Above  30  keV  we 
have  averaged  the  charge-exchange  cross  sections  of  Reference  73  for 
0*  + and  0*  + 0^  (reaction  4-10). 

0*  + N -►  0 + N*  and  NO*  » N ->-  NO  + N*:  Our  literature  search 

reveals  no  measurement  of  the  cross  section  for  either  of  those  reactions. 
For  ground-state  reactants  the  endothermicit ies  are  0.92  and  5.27  eV, 
respectively,  giving  threshold  speeds  of  4.9  x 105  and  1.0  * 106  cm/sec. 

NO*  + 0 ->  NO  + 0* : For  ground-state  reactants  this  reaction 

has  an  eiulothennicity  of  4.35  eV,  giving  it  a threshold  speed  of  9.0  * 10‘ 
cm/sec.  The  only  known  measurement  of  the  cross  section  (Reference  88)  is 
probably  strongly  affected  by  the  presence  of  excited  states  in  the  beams. 
Our  rough  estimate  is  that,  for  speeds  above  the  threshold,  the  cross 
section  equals  that  discussed  above  for  0?*  + O.  The  measurements  of 
Reference  88,  which  include  both  these  cross  sections  for  ion  energies 
from  70  to  6000  eV,  are  consistent  with  our  estimate  within  experimental 
uncertainty. 


Notes  to  Figures  4-1  to  4-8:  The  six  cross  sections  for  symmetric  charge 
exchange  (Figures  4-1  to  4-4)  and  the  eleven  for  asymmetric  charge  exchange 
(figures  4-5  to  4-8)  are  described  above  in  Section  4.  These  charge- 
exchange  cross  sections  arc  recommended  values  based  on  a literature  search, 
and  they  lean  heavily  on  experimental  rather  than  theoretical  cross  sections. 
They  arc  plotted  in  the  form  of  effective  momentum- transfer  cross  sections, 
q*(v)  or  q,.,(v)  (F:q,  2-43  or  8-31).  The  cross  sections  are  intended  to 

refer  to  ground-state  reactants  within  the  stated  accuracy;  comments  about 
the  states  of  the  reactants  in  the  experiments  have  been  included  above, 
for  symmetric  charge  exchange  the  recommended  ci*oss  sections  are  also  tabu- 
lated in  Table  4-2. 

Also  shown  on  figures  4-1  to  4-8  arc  momentum- transfer  ci'oss 
sections  for  elastic  scattering  qj,(v).  For  center-of-mass  energies  below 
20  eV  those  arc  calculated  from  F.q.  3-15  and  Table  3-2  or  arc  described  in 

O 

Section  3.3b.  For  high  energies  up  to  v - 10  cm/sec  the  cross  sections 
are  described  in  Section  3,3c.  In  Figures  4-5  to  4-7  we  have  interpolated 
a dotted  line  to  join  the  low  and  high  energy  portions  of  q^ (v ) . At  the 
highest  energies  q (v)  is  extrapolated  by  a straight  dotted  line, 

‘lD(v)  “ v~'\ 

The  momentum- transfer  cross  sections  for  inelastic  scattering 

q.  , (v)  arc  calculated  from  F.q.  0-8  and  Table  0-1.  Symbols  are  defined 
mel 

in  the  Glossary  in  Appendix  A. 

The  weight  of  a line  is  an  indication  of  accuracy.  A solid  line 
indicates  an  accuracy  within  +30  percent,  a dashed  line  lias  an  accuracy 
of  roughly  a factor  of  two,  and  a dotted  lino  is  for  order- of- magnitude 
accuracy.  There  are  two  kinds  of  dashed  lines  having  the  same  meaning; 
both  have  dashes  three  times  as  long  as  the  intervening  spaces.  Similarly 
the  dotted  lines  and  very  short  dashes,  which  are  only  as  long  as  the 


Table  4-2.  Total  cross  sections  for  symmetric  charge  exchange 
of  atmospheric  species.  The  values  are  best 
estimates  based  on  a review  of  the  literature. 


Relative 
speed  v, 
cm /sec 


Charge-exchange  cross  section  ox  (10  ^ cm2)  

ne+  + lie  N+  + N 0+  + 0 N;>+  + ^ °2+  + °2  N°f  + N0 


2(5) 

2.5(5) 

5(5) 

7(5) 


1(7) 

12.5 

20.5 

17.5 

2(7) 

10 

16 

14.4 

5(7) 

6.0 

11 .5 

10.8 

1(3) 

4.7 

8.5 

8.0 

1.5(3) 

3.2 

6.4 

6.3 

2(0) 

2.15 

4.5 

4.8 

3(3) 

0.93 

1.73 

2.65 

4(3) 

0.39 

0.71 

1.35 

5(8) 

0.17 

0.31 

0.7 

6(8) 

0.080 

8(3) 

0.017 

2(5)  means  2 * 10*  cm/sec.  The  number  of  digits  given  is  not  a 
measure  of  accuracy;  see  text  and  figures  4-1  to  4-4. 
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Momentum-transfer  cross  sections  for  elastic  scattering  q^v),  for  inelastic  scattering  q 
and  for  asymmetric  charge  exchange  q,(v)  of  0+  + 0„,  Of  + 0,  and  CL  + (L. 


ENERGY  OF  H ION 


RELATIVE  SPEED,  v (cm/sec) 

Momentum-transfer  cross  sections  for  elastic  scattering  qp(v),  for  inelastic  scattering  q 
and  for  asymmetric  charge  exchange  qx(v)  of  N+  + NO,  5+  + NO,  N0+  + N,  and  N0+  + 0. 


SECTION  5 

COEFFICIENTS  FOR  COLLISIONAL  MOMENTUM  AND  HEAT  TRANSFER 
FOR  TWO-FLUID  NUCLEAR-BURST  SIMULATIONS 


5.1  COUPLING  AND  HEAT-TRANSFER  COEFFICIENTS 

Sections  3 and  4 give  cross  sections  as  a function  of  the  relative 
speed  of  the  pair  of  particles.  We  now  integrate  these  over  a distribution 
of  speeds  to  obtain  averaged  quantities,  the  momentum-coupling  coefficient 
and  heat-transfer  coefficient,  as  a function  of  the  temperatures  and  macro- 
scopic velocities  of  the  two  gases.  These  arc  found  by  using  liqs.  2-33 
and  2-34  to  Integrate  the  momentum-transfer  cross  sections.  We  have  already 
done  these  integrals  analytically  for  two  cases:  in  Section  3.2  for  electron- 
ion  elastic  scattering  and  in  Section  5.3a  for  ion-neutral  elastic  scatter- 
ing by  the  polarization  force  alone.  We  now  integrate  numerically  for  the 
cases  of  electron-neutral  scattering,  ion-neutral  elastic  scattering,  and 
symmetric  charge  exchange. 

Tito  specific  cross  sections  treated  are  as  follows:  The  coeffi- 

* * + 
cients  s,  and  h,  due  to  symmetric  charge  exchange  (such  as  N + N -»■ 

4- 

N + N ) ire  based  on  the  cross  sections  in  Figures  4-1  to  4-4  (or  Table  4-2). 

(Asymmetric  charge  exchange,  such  as  n'*  + 0 ->  U+  + N in  Figure  4-5,  is 

omitted  here;  this  should  be  treated  as  a chemical  reaction,  ns  in 

Appendix  B.)  The  coefficients  s. ^ and  h . ^ for  elastic  scattering 

of  ions  and  neutrals  are  based  on  the  cross  sections  in  Figures  4-1 

to  4-8.  For  pairs  for  which  only  the  polarization-force  cross  section 

is  shown,  the  coefficients  s.,  and  h . , are  taken  from  Table  3-3,  while 

jk  jk 

for  the  other  five  pairs  (in  Figures  4-5  to  4-7)  they  are  found  by  numerical 
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for  elect  roil -neutral  scatter- 


integration.  The  coefficients  s.  and  h. 

j e j e 

ing  are  based  on  the  cross  sections  of  Section  3.1.  The  coefficients  s 

Kv 

and  h,  for  electron-ion  elastic  scattering  are  given  by  l;qs.  3-7  to  3-9 


and  3-5  with  2.n(2/0;|))  = 12.  Since  the  electron-ion  coupling  coefficient 


of  Hq.  3-7  is  essentially  independent  of  ion  mass,  it  equals  the  combined 


coupling  coefficient  Sje  of  L-q.  2-50. 


These  coupling  coefficients  and  heat-transfer  coefficients  arc 
found  by  integrating  Eqs.  2-33  and  2-34  analytically  or  numerically . The 
numerical  integrations  used  Simpson's  rule  with  a step  size  of  = 0.05 
over  the  interval  0 < £ < (V/U+4) ; the  resulting  integrals  are  accurate  to 
three  digits  or  more.  If  a cross  section  is  known  only  over  a range 
Va  < v < v , we  calculate  the  coupling  and  heat-transfer  coefficients  only 
for  temperatures  such  that 


1 

2 


v 


2 

Z 


nij  m2  2 u 


(5-1) 


The  coupling  coefficients  are  first  calculated  for  zero  relative  gas 

-y  -y 

velocity  = 0)  as  a function  of  the  weighted  mean  temperature 

(m^kTj+n^kTp/ (nij+n^)  (bq.  2-21)  from  0.01  to  500  eV.  figure  5-1  shows  the 
resulting  ion-neutral  and  electron- ion  coupling  coefficients.  The  ex- 
tension of  these  results  to  small  nonzero  relative  gas  velocity  is  discussed 
later  in  this  subsection.  Because  of  crowding  figure  5-1  docs  not  show  the 
(constant)  coupling  coefficients  for  0 + 0,  N + NO,  or  0 + NO;  these 

are  given  in  Table  3-3. 


The  accuracy  of  a coupling  coefficient  is  indicated  roughly  by  the 
weight  of  the  line.  The  solid  (most  accurate),  dashed  and  dotted  portions 
of  the  cross-section  curves  (figures  4-1  to  4-8)  have  been  repeated  in 
figure  5-1.  The  remarks  about  the  throe  weights  of  curves  made  in  Section 
4 apply  here. 
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gend 


Figure  5-2  shows  the  electron-neutral  coupling  coeffi clouts  and 
an  extension  from  Figure  5-1  of  the  electron-ion  coefficient,  which  is 
independent  of  ion  mass.  Figures  5-3  and  5-4  show  the  heat-transfer  coeffi- 
cients for  zero  gas  velocities.  These  are  related  to  the  coupling  coeffi- 
cients of  Figures  5-1  and  5-2  by  l-q.  2- 34b.  Because  of  crowding  Figure  5-3 
does  not  show  the  boat-transfer  coefficient  for  0 + N,  which  is  slightly 

greater  than  that  for  N+  + 0 (see  Figure  5-1  or  Table  3-3),  or  that 
for  N0+  + 0,  which  is  given  in  Table  3-3. 


We  have  also  computed  all  these  coefficients  for  nonzero  gas 


velocities,  jv^  , 


-v„ 


>0.  Presenting  the  results  is  simplified  by  using  anew 
independent  variable.  We  lot  v2  denote  the  mean  square  relative  speed  of 
pairs  of  particles  of  species  1 and  2;  the  new  independent  variable  t is 
then  (m  m,/(m1+m7).|v2/3,  which  is  2/3  of  the  mean  kinetic  energy  of  pairs 
of  particles  in  their  contor-of-mass  frames: 


m., 


T = 


— k i . 
(iii,+n>2)  1 


(iii. +m2) 


1 V? 

l-  '1‘  + — , — ' ■ 

K 2 3 (iiyra^,) 


YV2 


(5-2) 


We  call  T the  "collision  temperature"  for  collisions  between  species  1 
and  2;  when  V(  = V,  it  is  simply  the  weighted  mean  temperature  ki'm 
(L;.q.  2-21).  For  electron-neutral  and  electron-ion  coll  isions  the  neutral 
and  ion  temperatures  arc  usually  small  compared  to  lc»  so  l.q.  3 

reduces  to 

T = kT  + — m V - V , 

e a c | e 2 1 

where  V,  is  either  the  velocity  of  the  neutrals  or  that  of  the  ions. 


Using  X as  one  of  the  indeponden’:  variables,  one  finds  that 
the  coupling  coefficients  and  heat- transfer  coefficients  are  nearly  nule 
pendent  of  jVj-V0|  when 

" , 2 . k'[,  + 

^ - I ' 

which  is  usual lv  satisfied  in  late-t imo  two-fluid  simulations. 


(5-4) 


aend 
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To  indicate  this  we  have  labeled  tiie  abscissas  of  Figures  5-1  to  5-4  as  c. 
For  relative  speeds  satisfying  inequality  5-4  the  calculated  coupling  and 
heat-transfer  coefficients  are  always  within  ±15  percent  of  the  values 
given  in  Figures  5-1  to  5-4.  Larger  relative  speeds  make  substantial 
changes  in  these  coefficients,  however,  so  these  Figures  and  Table  5-1 
should  be  used  only  for  speeds  that  satisfy  inequality  5-4. 


Table  5-1  gives  empirical  fits  to  the  momentum- coupling  coefficients 
of  Figures  5-1  and  5-2.  Since  these  are  valid  only  for  low  relative  gas 
velocity,  heat- transfer  coefficients  can  be  calculated  from  them  by  the 
relation  h = 3s/(.m  +m2)  (Lq.  2-34b). 


5.2 


HOW  TO  USE  THE  COUPLING  AND  HEAT-TRANSFER 
COEFFICIENTS  IN  TWO-FLUID  SIMULATIONS 


We  have  presented  a method  of  treating  the  momentum  and  energy 
transfers  due  to  clastic  scattering  and  symmetric  charge  exchange  m ion- 
ncutral , electron- neutral , and  electron -ion  encounters.  We  give  a brief 
summary  for  the  reader  who  is  primarily  interested  in  results. 


The  neutrals  are  composed  of  several  species  (j=l,2,...)  having 


o 


particle  masses  in.,  number  densities  n'5.  and  total  number  density 
^ J 

n " n]  ■ The  neutrals  are  characterised  In  a mean  vclocitv  V and  a 
J J n 


n 


single  temperature  T ; we  consider  the  paiv.'le  velocity  distribution 


of  each  neutral  spec < to  be  Max.oji  1-Bolt ’.mann  in  the  frame  moving  at 
tiie  mean  velocity  V ,.ie  ions  are  composed  of  singly  charged  species 


'.L=  1 , 2 , . . . ) having  number  densities  n and  total  number  densitv  n = 

, k i 


i.,  n 


kV 
cJ  by 


• he  i-  and  electron  •ucio.ity  distributions  are  similarly  cliaracter- 


i , V_  and  T . The  electron  number  density  is  given  bv 


•t.  '*  n . . a no  t 1;  a mean  elect  ron  velocity  cun  ••e  e ! :i  i,> . ;<af  e<i  by  Fq.  2-46, 


'e.i 


(5-5) 


Table  5- 1 a . Parameters  of  the  coupling  coefficients  due  to  electron- 

neutral  scattering  Sje  and  elastic  electron-ion  scattering 
s-je  as  given  by  s = Br*3.  The  collision  temperature  x 
(Eq.  5-3)  is  in  eV  and  s is  in  gm-cnvfysec. 


sje  - Bt 


Pair  of 
species 

B 

b 

Range  ot  x, 
in  eV 

e + He 

5.90(-35)a 

+0.564 

0.01<t<0.8 

5.52 ( -35) 

+0.267 

0.8<t<4 

8.0(-35) 

0. 

4<t<10 

e + N 

3. 09 ( -35) 

+1.258 

0.3<x<2 

5 . 24  < - 35) 

+0.498 

2<x<10 

e + 0 

5.0(  -35) 

+0.530 

0.03<x<i 

5.0( -35) 

+0.748 

1 <T  <10 

e N2 

1 .43 (-34) 

+0.846 

0.01<x<0.8 

1 .30( -34) 

+0.436 

0.8<t<24 

b.2( -34) 

0. 

24<t<100 

e + 02 

8.30( -35) 

+0.903 

0.01<t<0.3 

5.42(-35) 

+0.549 

0 . 3<r<l 8 

1 ,27( -34) 

+0.255 

18<t<100 

e + NO 

5. C( -36) 

0. 

0.01<t<0.03 

7.27 (-34) 

+1.420 

0.03<x<0.2 

8.72(  -35) 

+0.102 

0.2<r<l  .7 

6.70(-35) 

+0.597 

1 ,7<x<13 

1 . 79( -34 ) 

+0.214 

1 3<t< 1 00 

e + ion 

3. 18(-32)e 

-1.5 

Note  d 

5.90  (-35)  means  5.90  x 10'°°. 
10'5  n01/3<  x <30  eV,  where  n. 


is  in  electrons/cnr . The  lower 


limit  is  from  Debye  theory,  as  mentioned  in  Section  3.3b. 

We  used  ln( 2/ 0m)  = 12  in  Eqs.  3-5  and  3-7.  This  factor  is  tabulated 
in  Reference  8,  Table  5.1  (called  in A)  and  plotted  in  Reference  7, 
Figure  8-6. 


Table  5-lb.  Parameters  of  the  coupling  coefficients  due  to  ion-neutral 
elastic  scattering  s^  and  to  symmetric  charge  exchange 
s£  as  given  by  sjk  = Bx^  and  s£  = B*rc.  The  collision 


temperature  x (Eq.  5-2)  is  in  eV  and  sj^ 
gm-cnr/sec. 


are  in 


Sect  ion  2. 3 gives  the  two-fluid  momentum  and  energy  balance  equa- 
tions in  a form  that  includes  elastic  scattering  and  symmetric  charge 
exchange  (such  as  N+  + N **■  N + N+).  We  repeat  these  relations  here  and 
then  summarize  how  to  calculate  the  nine  coefficients  for  collisional 
momentum  and  heat  transfer  (Eqs.  5-12  to  5-19).  Other  symbols  are  de- 
fined in  the  Glossary  in  Appendix  A.  The  momentum  balance  equations  for 
the  neutral  fluid,  for  the  ion-electron  plasma,  and  for  the  electrons 
(Eqs.  2-51  to  2-53)  are 


3(p  V ) + + 

!L2!_  + V • (p  v V ) 

3t  n n n-' 


- VP  + p g - n n.s  . (V  -V.)-n  n s (V  -V  ),  (5-6) 
n n&  n i nil  n l n e ne v n e" 


3(P;V.) 


+ V • (p.V.V.  ) = - V( P.+P  ) 4 X B ♦ (n.+0  )g 
3t  ill  i eJ  v'  i • e'& 


0 = - VP  - cn 
e e 


* \V,»  <VV  * VW<VV  - 

/ v 4\ 

!+  el  ->  ■>  -y  -> 

\E  + 1 + p g + n n s (V  -V  ) + n.n  s.  (V.  -V  ) 

\ cl  c n e ne  n e i e iev  i e' 


. (5-8) 


The  internal  energy  balance  equations  for  the  neutral  luid,  for  the  ions 
and  for  the  electrons  (liq,  2-55  to  2-57)  are 

-■r-  n)  + V*  (p  I V ) = - P V*V  + n n.s  ,d(n)  4.  2 

3t  n n n n n n l ni  in  n l 

Ml  n s dv'u  V -V  + n n.h  . k (T.  -T  ) + n n h k(T  -T  ) , (5-9) 

nenene  ne  nim  in  n e nc  v e n 1 


-►  -*•  (ni  ■+  ->  1 2 

— i i-r~-  + V*  (p.l.V. ) = - P.V*V.  + n n.s  . (l-din')  V -V.  r 
3t  xii  i i n i m in  n 1 1 

. ( i i -►  -t-  2 

+ nineSic  ic  VVe  " “ Vihnik  ~ Vehiek  (W  ’ t5'10) 
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+ V-  n kT  V I = - l>  V • V 
\2  q e c)  e e 


+ n n s 

i-d(n])  v 4 

2 / 

i n.n  5, 

I -df^) 

V.-tf 

n e ne  \ 

ne  / n e 

i e te  \ 

ic  f 

l e 

2 


- n n h k(T 
n e ne 


+ n.n  h.  k(T.-T  ) 
i e le  v i e 


(5-11) 


These  six  balance  relations  should  include  additional  terms  to 
account  for  chemical  reactions.  (Asymmetric  charge  exchange,  such  as  re- 
actions 4-7  to  4-20,  is  classed  in  this  work  as  a chemical  reaction.  See 
Appendix  B and  Figs.  4-5  to  4-8.)  liq.  5-11  should  also  include  an  addi- 
tional term  to  account  both  for  the  heat  transfer  resulting  from  vibrational 
excitation  and  de-excitation  of  N,  by  electron  impact  and  for  the  electron 
energy  loss  resulting  from  electronic  excitation  and  ionization  of  neutrals 
and  ions  by  electron  impact. 


The  coupling  and  heat-transfer  coefficients  depend  on  the  fluid 


composition.  The  first  one,  s 


ni 


is  calculated  from  Eq.  2-48  as  the 


weighted  sum  of  the  coupling  coefficients  s.,  for  elastic  scattering 

j K 

of  pairs  of  species  consisting  of  a neutral  species  j and  an  ion  species 

k (N+  + N,  N+  + 0 N+  + N.)(  etc.)  plus  a sum  of  the  coupling  coefficients 

s*  due  to  symmetric  charge  exchange  of  ions  with  their  parent  atoms  or 
^ + + 

molecules  (N  + N,  NO  + NO,  etc.);  the  relation  is 


n n . s . 
n l ni 


EL 

.i  k 


o ♦ 
n.  n. 
) k 


S ...  + 
jk 


L 

k 


o + * 

nk  nk  sk 


(5-12) 


where  n?  is  the  number  density  of  neutral  species  j and  nj 
is  the  number  density  of  ion  species  k.  The  coupling  coefficients  s . 
and  s*  for  individual  pairs  of  species  can  be  found  from  Table  5-lb  as 
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a function  of  the  collision  temperature  c (Kq.  5-2)  for  that  pair. 


Table  5- lb  applies  when  the  ion-neutral  relative  speed 
inequality  5-4. 


V -V. 
n i 


satisfies 


The  ion-neutral  heat-transfer  coefficient  h . should  be  summed 

m 


at  the  same  time;  it  is  given  by  Kq.  2-59, 

Vihni  * S?  njnk  Ifrib)  + ? nknk  2m 


■?s* 
5 k 


(5-13) 


j k 


■j  "V 


(») 


for  the  fractional  coupling  coefficient  d^. ' one  may  use  the  approximate 
Lq.  2-66  rather  than  sum  with  the  definition,  Kq.  2-53a;  these  relations 
are 


n n.s  ,d ^ e V T]  n°n+  — — -?■  s..  +■ 
n i iu  ni  t-r*  y j k (in.+m.  ) jk 

J K j K 


V*  o+l 
V ' nvui-  7 


k k 2 °k 


0.4  n n.s  . 
n i ni 


(5-14) 


The  electron-neutral  coupling  coefficient  s is  calculated 

1 no 

from  Kq.  2-49  as  a sum  over  all  neutral  species: 


V'  o 

n n s = n >,  n s. 
n e ne  e " j ie 


(5-15) 


The  individual  coefficients  s.  are  found  from  Table  5-la  as  a function 

.1° 

of  the  electron-neutral  collision  temperature  T,  which  is  given  by  Kq.5-3. 
The  clcctron-neutral  heat-transfer  coefficient  h is  given  by  Kq.  2-60, 


n n h ^ n 21  n°  3s.  /m, 
n e no  e l ie 

J 


(5-16) 


which  can  be  approximated  by  Kq.  2-62.  The  coefficient  d^J^  is  defined 
by  Kq.  2-54,  but  one  can  usually  approximate  by  using  the  mass  of  the  N., 
molecule  (Kq.  2-63);  these  relations  are 
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, in  i n u v 

n n s u _ n n • r s . 

n c ne  ne  ~ e . j (m.+m  ) ie 
J J oJ 


2.0  ’ 10  ^ n n s 

n o ne 


(5-17) 


The  clectron-ion  coupling  coefficient  s . ^ , which  is  independent 

of  the  composition  cf  the  ion  fluid,  is  given  by  the  bottom  line  of  Table 

5-la  or  by  ! lq s . 5-7,  5-5,  and  5-9.  It  is  a function  of  the  electron-ion 

collision  temperature  r,  which  is  given  by  Eq.  5-5.  The  heat-transfer 

coefficient  h.  is  given  by  Kq.  2-61  and  normally  can  be  approximated 
1 e + 

by  using  the  mass  of  the  0 ion  (Eq.  2-64);  these  relations  arc 


n.n  h.  i 3n  i.  V iC/m, 
j.  e re  e ie  ^ k k 


5n ,n  s . /m (0  ) 

1 c ie  v ' 


(5-18) 


The  coefficient  df^  involves  the  same  sum  and  is  given  by 
re  * 


, (i) 

n . n s . d . ' = n m s 


l e ie  ie 


V,  n,  /m. 
c c ie  , k k 


3.4  x 1()  n.n  s. 

i e l e 


(5-19) 


Our  coupling  coefficients  may  be  more  familiar  if  they  are  in- 
terpreted as  collision  frequencies  and  related  to  the  electrical  resistivity 
The  average  momentum- loss  collision  frequency  for  electron- ion  collisions  is 


ei 


n . 

l 

m 

e 


s . 

ie 


(5-20) 


and  that  for  collisions  of  electrons  with  neutrals  of  species  j is 


0 

n , 

V • = — S. 

ej  me  je 


(5-21) 


These  quantities  arc  related  to  the  electrical  resistivity  n (in  cm^/sec)  by 


c m 

= -T-(vci  *?v  A 

e % \ J / 

= 4 /s.  ♦ ± Z «°s.  ) 
e*  \ ie  ne  j •> 


(5-22) 


The  electron  momentum  balance  equation  (r.q.  2-53  or  5-8)  can  be 

expressed  in  terms  of  the  resistivity  by  using  l'.qs.  5-22,  5-5,  and  2-49. 

If  one  can  neglect  the  magnetic  field,  the  electron  pressure  gradient, 

gravity,  and  the  ion-neutral  slip  velocity  V -V.,  the  resulting  relation 

-►  ->  " * 

reduces  to  a form  of  Ohm's  law,  E * hJ/c.  More  generally,  the  resulting 

•> 

relation  can  be  substituted  into  Faraday's  law,  DU/D t = - cvxH,  to  give 

5-5-  = V * j V.xn  - nj  - SHI  + - £ n s (V  -V.)l  . (5-25) 

^ i cir  ene  e n nev  n l'j  v 

Mill)  simulations  usually  ignore  the  last  tin-ee  terms,  although  the  validity  of 
this  approximation  is  uncertain. 

Simplified  forms  of  the  momentum-coupling  coefficients,  heat- 
transfer  coefficients,  and  electrical  resistivity  n are  useful,  for  the 
temperature  ranges  and  species  mixtures  encountered  after  the  first  second 
of  a nuclear  bursts,  the  following  forms  are  within  a factor  of  two  of  re- 
sults found  by  using  Eqs.  5-12  to  5-22  with  Table  5-1  or  Figs.  5-1  to  5-4: 

n n.s  . = 1.2*  10  { (N*0)  «■  1.6  [N,+0,+N0] } f N++0%N0+] 

n i m 2 2 
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(5  - 24)  * * 
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(5-25)  * * * 
(5-2t>)  ' * * 
(5-27) * * * 
(5-28)  * * * 
(5-2i>)  * * * 
(5-30)  * * * 


where  kT  , kT.,  and  kT  are  in  eV,  n is  in  cm"/ see,  and  the  other 
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quantities  are  in  cgs  units.  In  liqs.  5-28  to  5-30  we  used  £n(2/0  ) = 12; 

Hi 

the  weak  dependence  of  this  factor  on  n^  and  T0  is  given  in  Reference  8 


Table  5.1  (where  £nA=f.n(2/0  ) ) or  in  Reference  7,  figure  8-0. 


Using  the  simplified  form  of  the  momentum-'  >upling  coefficients, 
the  ion-neutral  collision  frequency  v.  , the  electron-neutral  collision 
frequency  v , and  the  electron- ion  collision  frequency  v^.  are  given  by 

/p.  = 1 .2x  |0~-l2{  |N+0]  + ! .f,[N  +0)+N0]  } [N++O+«N0+]/p. 

i £.  m.  1 

* 4. 8x|  O' '^(kT  +kT.  )°"'/5U.2|N|  |N*  ] + |0|[0+|  a 0. 84  | NO  | |N0+  | }/p.  , 
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n e ne  e e n 


\>  . = n.n  s.  /p  = 3.5x10  (kT  ) '"'’n., 
ci  i e ie  e e i 


(5-31 1 -*  * 

where  kT  . kT. , and  kT  are  in  oV  and  the  other  quantities  are  in  cgs  units, 
n l e 1 
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APPENDIX  A 
GLOSSARY  OF  SYMBOLS 


(W.S) 


F(W) 

5 


Bohr  radius,  0.5292  * 10"^  cm. 

Constant  in  electron-ion  (Coulomb)  scattering  cross  section, 
0.06  x 1017  l'n(2/Om)  cm6/sec4  (Eq.  3-5). 

Magnetic  field  strength  (gauss). 

Speed  of  light,  2.998  * ]0^  cm/sec. 

Fractional  coupling  coefficient:  d^l  fraction  of  the  ion- 
neutral  frictional  heating  that  goes  into  internal  energy 
of  neutrals,  Eq.  2-53;  (1^1  d^l  fractional  coupling  coeffi- 
cients for  electron-neutral  and  electron-ion  collisions, 

Eq.  2-54. 

Proton  charge,  4.803  * 10" ^ statcoul. 

Base  of  natural  logarithms. 

Electric  field  strength  (statvol t/cm) . 

Maxwel 1-Bol tzmann  velocity  distribution  function,  Eq.  2-1; 
f^(v-j),  abbreviated  form  defined  by  Eq.  2-16;  f2(v),  distri- 
bution of  relative  velocities  of  pairs  of  particles  of 
different  species,  Eq.  2-17  (sec^/cm^). 

Dimensionless  momentum- transfer  cross  section  for  elastic 
atom-atom  scattering,  q^/\  (Eq.  3-25). 

2 

Acceleration  due  to  gravity  (-980  cm/sec  ). 


1 


8-1 


h(T(||,V) 


k(Tnl,v) 


K0(T) 


«.  (T  ,V) 
x'  m ' 


Planck's  constant,  6.626  - lO-*^  erg-sec. 

Rate  coefficient  for  heat  transfer  due  to  scattering,  Eq.  2-30 
(see  also  Lq.  2-40);  h^,  heat- transfer  coefficient  for  elastic 
scattering  of  neutrals  of  species  j and  ions  of  species  k; 
hjfi,  hj.^  (or  h-e  in  Section  3.2),  heat-transfer  coefficient  for 
scattering  of  electrons  and  neutrals  of  species  j or  ions  of 
species  k;  h.^,  hj.e,  heat- transfer  coefficient  for  scatter- 
ing of  electrons  and  neutrals  of  species  j or  ions  of 
species  k;  h£,  heat- transfer  coefficient,  due  to  symmetric 
charge  exchange  of  ions  of  species  k with  their  parent 
atoms  or  molecules,  Section  2.2b;  h^.,  hne>  h • e , combined 
heat- transfer  coefficient  for  ion-neutral,  electron-neutral, 
or  electron-ion  collisions.  Section  2.3b  (cm  /sec). 

Specific  internal  enerq.y  of  a gas,  Section  2.1a.  This  is 
taken  to  include  the  translational  kinetic  energy  of  particle 
motion  in  the  frame  moving  with  the  mean  velocity  plus 
the  energy  of  rotational  excitation  of  molecules.  Vibrational 
energy  and  ionization  energy  are  excluded,  (erg/gm) 

Index  of  neutral  species. 

Current  density  (abamp/cm4) . 

Boltzmann  constant,  1.381  * 10  ^ erg/°K  or  1.602  * 10*^  erg/eV. 
Index  of  ion  species. 


Reaction  rate  coefficient:  ks,  rate  coefficient  for  scatter- 
ing collisions,  Eq.  2-28,  used  only  for  inelastic  collisions; 

k , rate  coefficient  for  asymmetric  charge  exchange,  Eq.  B-8 

* 3 

(see  also  Eq.  B-14)  (cm  /sec). 

Mobility  of  ions  in  a weak  electric  field,  reduced  to  standard 

p 

density  (cm  /statvol t-sec) . 

Rate  coefficient  for  specific  internal  energy  change  of  a 
reactant  due  to  asymmetric  charge  exchange,  Eq.  B-10  (see 
also  Eq.  B-20)  (cm-sec). 

A- 2 


P (T  , V) 
rx  111 


S<VV> 


Particle  mass;  in  , electron  mass;  m, , mass  of  a molecule  of 
neutral  species  j;  m^  or  nr,  mass  of  an  ion  of  species 
k (gm) . 

Number  density  of  particles;  n , n.,  n , number  density  of 

^ " o Q 

electrons,  ions,  or  neutral  molecules;  n(X.)  or  n.,  number 

J J 

density  of  neutral  species  j;  n(X^)  or  n*,  number  density 
of  ions  of  species  k (cm‘2). 

date  coefficient  for  velocity  change  of  a reactant  due  to 
asymmetric  charge  exchange,  Eq.  B-9  (see  also  Eq.  B- 1 9 ) 

(cm-sec) . 

fluid  pressure;  Pe,  P. , Pn,  pressure  of  electrons,  ions,  or 
neutrals  (dyne/cm^). 

Momentum-transfer  cross  section  of  a pair  of  particles, 

Eq,  2-23;  q^(v),  momentum-transfer  cross  section  for  elastic 
scattering;  M^ne](v),  momentum- transfer  cross  section  for 
inelastic  scattering,  Eg.  C-3;  q*(v),  effective  momentum- 
transfer  cross  section  due  to  symmetric  charge  exchange, 

2 o (v)  (Eq.  2-43);  qT(v^ ) , effective  momentum- transfer  cross 
section  due  to  asymmetric  charge  exchange,  ((m,+m?)/m, ) a (v,) 
(Eq.  B- 31 ) (cm2). 

Separation  of  two  particles;  r , separation  at  which  their 
potential  energy  is  a minimum  (cm). 

Coupling  coefficient  due  to  scattering,  Eq.  2-29  (see  also 
Eq.  2-36);  s^,  coupling  coefficient  for  elastic  scattering 
of  neutrals  of  spec'es  j and  ions  of  species  k;  s-  , 

J ^ 

coupling  coefficient  for  scattering  of  electrons  and  neutrals 
of  species  j;  s^,  coupling  coefficient  for  elastic  scattering 
of  electrons  and  ions  of  species  k;  s*,  coupling  coefficient 
due  t:o  symmetric  charge  exchange  of  ions  of  species  k with 
their  parent  atoms  or  molecules.  Section  2.2b;  s ^ , $ne,  s^e, 
combined  coupling  coefficient  for  ion-neutral,  electron  neutral , 
or  electron-ion  collisions,  Eqs.  2-48  to  2-50  (gni-cnr /sec) . 


Stopping  cross  section,  -n^_^dC-j/dx  (r.q.  C - 1 ) ; S^nei(v|)> 

stopping  cross  section  duo  to  inelastic  scattering  from 

2 

atomic  electrons  (erg-cm  /atom). 

Time  (sec). 

Temperature;  T0,  Tj,  Tft,  temperature  of  electrons,  ions,  or 
neutrals  (°K  or  oV ) . 


Temperature  characterizing  the  distribution  of  relative 
velocities  of  pairs  of  particles,  (ii^Tp  + nuTj )/(m^+m^) 

(Eq.  2-21);  T^,  temperature  characterizing  the  distribution 
of  relative  velocities  of  pairs  consisting  of  a neutral  atom 
or  molecule  of  species  j and  and  ion  of  species  k (“K  or  eV). 


Thermal  speed  parameter,  (2kT^/m^  + ZkTg/m,,)^,  character- 
izing the  distribution  of  relative  velocities  of  pairs  of 
particles  (Eq.  2-32)  (cm/sec). 


Initial  relative  velocity  of  two  particles,  v(-Vr,  (Eq.  2-6); 
v',  final  relative  velocity  of  two  particles;  v^.,  typical 
relative  speed,  Eq.  3-4  (cm/ sec). 


Initial  velocity  of  first  particle  in  the  laboratory  frame; 
v{ , final  velocity  of  first  particle  (cm/sec). 

Velocity  of  center  of  mass  of  two  particles  (cm/sec). 

Relative  macroscopic  velocity  of  two  gases,  U, -V^  (EQ-  2-20) 
( cm/ sec ) . 


Macroscopic  velocity  of  first  gas  in  the  laboratory  frame; 
\?e,  , l?n,  macroscopic  velocity  of  electrons,  ions,  or 

neutrals  (cm/ sec). 

Dimensionless  kinetic  energy  of  a pair  of  particles  in  their 
center- of -mass  frame,  Eq.  3-24. 


Atomic  number. 


Anqle-averatjed  dipole  polarizability  (civ). 

Dimensionless  parameter  measuring  the  relative  strength  of 
the  r*k  and  r"^  potential  energies.  Eg.  3-18. 

Dirat  delta  function. 

Change  in  velocity  of  first  particle  during  a collision, 
v-j  *-v-|  (cm/sec). 

Rate  of  change  of  mass  density  due  to  asymmetric  charge 
exchange,  f.q.  B-2  (gin/ cm  -sec). 

Rate  of  change  of  internal  energy  density  due  to  scattering 

or  symmetric  charge  exchange,  Eg.  2-38;  /.(pi ) /At,  rate  due 

3 * 

to  asymmetric  charge  exchange  (erg/cm  -sec). 

Rate  of  change  of  total  energy  density  due  to  scattering  or 
symmetric  charge  exchange,  Eg.  2-3;  a|]>  p v^j  /At,  rate  due 
to  asymmetric  charge  exchange.  Eg.  8-4  (erg/ cm3 -sec) . 


Rate  of  change  of  momentum  density  due  to  scattering  or 

symmetric  charge  exchange,  Eg.  2-2;  A(p’v)  /At,  rate  due  to 

x 3 

asymmetric  charge  exchange.  Eg.  B-3  (dyne/cm  ). 

Loss  (or  gain,  if  r<0)  of  kinetic  energy  by  a pair  of 
particles  in  a collision;  e^,  loss  of  kinetic  energy  in  a 
collision  in  which  the  final  states  of  the  particles  are 
specified  by  v (erg). 

Depth  of  potential  energy  well  of  a pair  of  particles 
(erg  or  eV)  (used  in  Section  3.3b). 

Scaling  factor  for  energy,  Eq.  3-23  (erg  or  kev). 

Electrical  resistivity,  (c  /e  )(s ie+ne"  jnjSj„)  Eg. 3- 22) 

n-.  resistivity  transverse  to  a strong  magnetic  fio'd, 
(4/3)(2mii  l^Vc"  t’ii(2/  ' )/(kl)3^3  (enf/sec,  emu). 


Zenithal  scattering  angle  in  the  center-of-mass  frame, 
defined  by  diayram  in  Section  2.1b;  " , minimum  scattering 
angle,  Cq.  3-2. 

Scaling  factor  for  length,  Cy.  3-22  (cm). 

Debye  shielding  distance,  (kT  /4-ie^n  (cm). 

Reduced  mass  of  a pair  of  particles,  m-m„/(m^  tm,,) 

(Cq.  2-4)  (ym). 

Average  momentum- loss  collision  frequency  for  electron-ion 
collisions,  Cq.  5-5;  average  mumentum-loss  collision 
frequency  for  collisions  of  electrons  with  neutrals  of 
species  j,  Cq,  5-6  (sec  ). 

Dimensionless  relative  velocity  of  two  particles,  Cq.  2-31. 

Mass  density;  , , , , , mass  density  of  electrons,  ions, 

(•  u 

or  neutrals  (gm/cm). 

Total  cross  section;  o t(v),  total  scattering  cross  section, 

Cq.  2-22,  used  only  for  inelastic  collisions;  o (v),  total 

2 x 

cross  section  for  charge  exchange  (cm  ). 

Differential  scattering  cross  section  in  the  center-of-mass 

O 

frame  (Reference  i)  (cm/sterad ) . 

Collision  temperature  for  collisions  between  two  species, 
Cqs.  5-2  and  5-3  (erg  or  eV). 

Azimuthal  scattering  angle  in  the  center-of-mass  frame, 
measured  around  the  v axis. 

Potential  energy  of  a pair  of  particles  (erg). 

Solid  angle  (storad). 
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APPENDIX  B 


MOMENTUM  AND  ENERGY  1RANSEER  BY 
CIIARGL-EXCIIANGC  REACT  IONS 


(a)  Statement  of  the  Problem 


Asymmetric  charge -exchange  react  ions  require  a more  general  treat- 
ment than  that  given  in  Section  2...’tt>)  for  symmetric  ones.  Ke  consider 
reactions  of  the  form 


a*  « is  • a i n* 


i ,1  ; ■’  I'5-1) 

we  will  sometimes  refer  to  the  species  l>>  the  numbers  given  under ncath  them, 
The  quantity  >-  again  denotes  the  loss  (or  gain)  of  kinetic  energy  of  the 
pair  of  particles,  which  for  endothermic  reactions  equals  the  ionization 
energy  ot  species  11  less  that  of  A fif  A and  1>+  are  in  their  ground 
states).  We  assume  that  this  is  negligible,  tluit  is,  |>|  < --  -pv“. 


We  further  assume  that  charge  exchange  results  only  in  forward 
scattering  fo  I);  this  approximation  is  good  at  energies  above  0.1  oY. 
While  it  would  be  easy  enough  to  account  for  non-forward  scattering,  we 
have  found  no  theoretical  or  experimental  studies  of  the  differential 
cross  section  for  charge  exchange  at  energies  below  0.1  eY. 


The  total  cross  section  for  charge  exchange  is  denoted  by  •»  (v) 


(analogous  to  o (v) , Lq.  2-22).  The  rates  of  change  due  to  charge  exchange 
of  the  densities  of  mass,  momentum,  and  total  energy  of  species  A*  are 


, ;j  ' 

- / m,  fj(Vj)f  ,lv  Jvox(v)d'\-jd->v, 


-=  • 7-  x- lx  t t .-wToi 


ssftia&iri iaesw 


A(*VVx  f -r  ■>  ■ r 3 s 

At  = • Jmlvl<'l(vl)t2(v2)v0x(v)d  Vld  v2  * 
A('2  Vi  )x  f 1 2..  , .3  .3 

At = J 2 Vl  fl(vl)f2(v2)vox(v)d  Vjd 

For  the  product  species  A the  corresponding  rates  are 


v2  • ( 


Ac.'  Ap. 

3x  lx 


A^3^3^  f ■>  •+  3 3 

— 1*~  = 1V1  f1(v1)f2(v2)vax(v)d  vxd  v2  , (B-6) 

— 2--k~~  "/  1^!  v{2f1(v1)f2Cv2)vax(v)d3vid3v2  • (r-7) 


(b)  Evaluation  of  Integrals 

The  procedure  for  evaluating  these  integrals  is  similar  to  that 
described  in  more  detail  in  Section  2.1(c).  We  change  integration 
variables  from  v^,  v2  to  v^,  v and  integrate  over  d v^  in  liqs.  B-2 
to  B-4.  The  integrals  over  d°v  are  then  expressed  in  terms  of  the 
following  three  rate  coefficients.  The  rate  coefficient  for  the  charge 
exchange  reaction  is  defined  analogously  to  Uq.  2-28: 


W>  EirA<;>Vv|v>*3 


Wo  also  define  two  related  coefficients  by 


w>  ? ; t<.  iJ  i2(oox(,)vc;-v,d3v , ;k- 


fat»»wsa»!«,-*: 
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Note  that  k is  in  cm  /sec,  while  p and  £ are  in  cm- see.  We  will 

show  that  p and  can  be  interpreted  as  the  rate  coefficient  for 

velocity  change  of  a reactant  and  that  for  specific  internal  energ'  change. 

We  do  the  angular  parts  of  these  integrals  as  in  Section  2.1(c). 

-V 

The  integral  in  liq.  11-9  is  parallel  to  the  V direction.  We  again  use 
*►  -*■  ->■ 

the  abbreviation  V = V,  - V,  and  eliminate  the  mean  temperature  T|n  m 
favor  of  the  related  thermal  speed  parameter  (liq.  2- .12) 


U.  P 

N "i 


2kT, 


The  result  for  k (T  ,V)  has  a form  identical  to  that  of  liq.  2-35  for 
xv  m 

Wi-- 

'i“i  fo  (U<K 2 fo-K-v/”)2  - c-«>v/u>2]  «H 

/iT  v J x L 


kx^Tm’V-* 
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tx 
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ox(UOf,Vf'  d£ 


if  T >0  and  V > 0, 
m 


if  V = 0 , 


Vox(V) 


if  T = 0 

nt 


(B- 11) 


The  results  for  the  other  rate  coefficients  are 
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T > 0 and  V > 0 , 

in 


2 


if  T >0  and  V = 0 , 

m 
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+ 2 -Am  trV  ^(V) 


if  T =0  and  V > 0 
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(c)  Results  for  Reactant  Species 

'Hie  time  rates  of  change  due  to  charge  exchange  of  the  densities 
of  mass,  momentum,  and  total  energy  of  the  first  gas  are 


Ap 


lx 


= - n.n.m.k  (T  ,V) 
At  1 2 1 m ' 


(B-14) 


A(PiVi)x 

At 


' nln2[kxClm’V)Vl  + Px(Tm’V)kTlV]  ’ (B'15) 
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. (B- 16) 


We  again  split  the  total  energy  into  macroscopic  kinetic  energy 
and  internal  energy  as  in  hq.  2-38.  The  macroscopic  kinetic  energy  is  now 
found  by 


4 ‘■l\2)x  i 


At 


At 


AWx  - _lv2Aplx 


At 


1 2 1 At 
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so  the  rate  of  change  due  to  charge  exchange  of  the  internal  energy  density 
of  the  first  species  is 

mp^Px  h (*V2  1 

— £r~  = - v2[f  kTikx(Tn.’v)  + WV>J ■ (B-18) 
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Hie  quantity  Px t 1 m * v ) wight  be  considered  to  be  the  rate  coeffi- 
cient for  velocity  change  of  a reactant  because 


1 

X 
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^iVx 

Ap, 

V lx 

At 

°1 

At 

1 At 

k'V 


= - n^p  (T  , V) 
2‘  x m ’ 


Ivv 


(B-19) 


This  rate  coefficient  PxCi’m>V)  can  be  positive  or  negative.  If  °x(v) 

is  inversely  proportional  to  v (say  o (v)  = K /v) , then  p (T  ,V)  and 

x o x ra  ~ .+. 

AV^/At  are  2ero  (see  definitions  B-9  and  2-17  of  PX(T„>V)  and  f2(v)); 

in  this  case  every  particle  of  the  first  species  has  the  same  probability 

per  unit  time  of  charge  exchange  (n^K^).  In  general,  PxO'm>V)  is  positive 
if  vax (v)  is  an  increasing  function  of  v (as  is  usual  for  symmetric 
charge  exchange,  such  as  Figures  4-1  to  4-4)  and  is  negative  if  vo^Cv) 
is  a decreasing  function. 


Similarly  we  note  that 

AI,..  , |A(p,I,)„  Ap 


lx 

At 


1 

PWx 

P1 

At 

lx 
At  . 
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--  % (T  ,V) 
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(B-20) 


1 


where  we  liave  taken  1^  = jkT^/m^.  One  could  call  &x(Tm>V)  a rate 
coefficient  for  specific  internal  energy  change  of  a reactant.  If  0 (v) 
is  inversely  proportional  to  v,  then  £x(i'm,V)  and  AI^/At  are  iero 
(see  definition  B-10). 


B-6 


(d)  Results  for  Product  Species 

Hie  relations  for  product  species  arc  quite  similar  to  those  for 

reactant  species  because  of  our  assumptions  mentioned  earlier  that  charge 

exchange  results  only  in  direct  forward  scattering  (0  <<  1)  and  that  the 

12 

energy  defect  is  negligible  (|ej  « 2"  yv  ).  In  this  case  the  rates  for  the 
product  that  has  particle  mass  simplify  to 


(B-21) 


A^'Vx 

At 


At 


v/l  _ a(?  Vi2),  , _ L 

At  At  (mj+m2)  m^  At 


(B-22) 


(B-23) 


The  last  term  of  bq.  B-23  represents  a part  of  the  energy  loss  (or  gain) 
occurring  in  each  reaction  (compare  liqs.  2-37  and  2-40);  it  need  be  included 
only  if  strict  energy  conservation  is  demanded. 

IVe  again  split  the  total  energy  P5V^2  into  macroscopic  kinetic 
energy  and  internal  energy  as  in  Bq.  2-58.  The  macroscopic  kinetic  energy  is 


A(  2 P3V3~)x  _ 1_  a[(P3v53  /p3]x 


A^3V5^x 


-+  1 2 Ap3x 

V . i V — 

3 2 3 At 
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for  the  internal  energy  of  this  product  species  the  rate  of  change  due  to 
charge  exchange  is 


. 4 P3V32)x  _4  l(p3V32), 


and  after  substitution  of  Eqs.  B-23,  11-24,  2-38,  and  B-17  this  becomes 


A(p  I,  A(p  I ) ^ A p V Ap 

At  At  *■  3 V At  2 1 1 3 •'  At 
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(m1+m9)  At 


(B—  25) 


The  total  power  per  unit  volume  transferred  by  charge  exchange 
from  the  total  energy  of  both  reacting  gases  to  ionization  energy  is 
(compare  liq.  2-41) 


4‘v/)*  4v/) 


jrz 3 n n.Ek  (T  ,V) 

At  1 2 x m ' 


(B-26) 


One  hydrocode  that  uses  a finite-difference  scheme  (the  MICE 
code)  uses  a special  technique  when  two  fluid  elements  having  different 
velocities  are  to  be  combined:  it  conserves  kinetic  energy,  as  well  as 

mass,  total  energy,  and  the  direction  of  momentum,  and  sacrifices  con- 
servation of  the  magnitude  of  momentum.  Use  of  this  technique  modifies 
Eqs.  B-24,  B-25,  and  B-22. 


(e)  A Simpler  Approximation 

We  present  a simplified  treatment  that  uses  only  the  ordinary 
rate  coefficient  for  a charge-exchange  reaction  instead  of  three  rate 
coefficients. 


Much  simplification  results  if  one  assumes  that  the  momentum  and 
internal  energy  densities  of  reactant  species  can  bo  approximated  by 


and 


A(plVx  v Aplx 
At  ' V1  At 


&frlVx 
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3 kFl  Aplx 


2 m. 
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(B-27) 


(B-28) 


The  neglect  of  terms  p^AV^/At  and  p^AI^/At.  in  these  relations  is 

equivalent  to  ignoring  the  rate  coefficients  p (T  ,V)  and  £ (T  ,V) , 

x ni  x in 

respectively  (Eqs.  B-19  and  B-20).  This  fact  enables  one  to  simplify  the 
other  relations  in  the  preceding  two  subsections. 


This  approximation  is  also  equivalent  to  assuming  that  the 
particles  that  charge-exchange  have  the  same  mean  velocity  and  specific 
internal  energy  as  the  rest  of  their  species.  The  direction  or  sign  of 
the  neglected  terms  can  be  scon  from  subsection  (c) : for  example,  if 
cx(v)v  is  an  increasing  function  of  v,  then  p^AV^/At  is  in  the 
direction  of  V?  - V^. 

(f)  Comparing  the  Effect  of  Charge  Exchange  with  that  of  Scattering 

finally,  we  ask  how  to  compare  the  effects  of  asymmetric  charge 
exchange  with  those  due  to  scattering  and  symmetric  charge  exchange,  first 
note  that,  the  former  causes  transfer  of  mass  between  ions  and  neutrals 
while  the  latter  two  do  not.  Then  for  comparing  momentum  transfer, 
consider  the  following  experiment. 

A beam  of  ions  of  one  species  is  incident  with  speed  on  a 

cold  stationary  gas  of  a second  species.  While  the  first  gas  traverses 
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4 


B-9 


ig 


iikaiiWHtti 


»«WvWw,  
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a distance  Ax  = VjAt,  charge  exchange  changes  Its  momentum  density  at 
a distance  rate 


A(Plvi) 

Ax  = ' "lmlVln2°x(Vl)  : 


(B-29) 


this  can  be  obtained  from  t:qs.  U- IS  and  U-llc  with  Tj  = 0,  - 0,  and 

V,  = 0.  for  comparison  scattering  off  the  second  gas  changes  this  momentum 
density  at  a rate  given  by  F.q.  2-36,  which  in  this  case  can  be  rewritten  as 

(«y»2)  A(PjVp 


n(vj)  - - Ax 


(B-30) 


This  could  be  considered  a definition  of  the  momentum-transfer  cross 
section  q, 

We  have  extended  this  relation  to  define  an  effective  momentum- 
transfer  cross  section  for  charge  exchange  q,(,(V^).  Using  Uq.  B-29  one 
finds  that 


VV  - 


(m1+m2) 


°x<V 


(B-31) 


This  definition  is  consistent  with  that  of  the  momentum- transfer  cross 
section  for  symmetric  charge  exchange  q*(v)  (liq.2-43).  The  relative 
importance  of  scattering  and  charge  exchange  in  transferring  momentum  is 
seen  by  comparing  q^,  with  q()  and  q^  ^ (see  Figures  4-5  to  4-8). 
Note  that  this  comparison  applies  only  to  the  momentum  change  of  the  ions 
because  definition  B-31  is  asymmetric  between  the  two  species.  uso  the 
comparison  applies  strictly  only  in  this  experiment,  in  which  Tj  = 0, 

T,  = 0,  and  V = 0. 
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APPENDIX  C 

INELASTIC  I ON -NEUTRAL  COLLISIONS  ABOVE  10  keV 


In  collisions  of  atoms  or  ions,  inelastic  scattering  involving  ex- 
citation  and  ionization  is  important  at  relative  speeds  above  5 * 10  cm/sec 
(centcr-of-mass  kinetic  energy  above  10  keV  tor  N ions  on  N) . After  a high- 
altitude  nuclear  explosion  such  collisions  are  usually  important  only  for 
the  first  second. 


We  use  measurements  of  the  energy  loss  rate  dU/dx 
for  beams  of  N and  M ioas  in  cold  N7  or  air.  The  measurements  of 
Itcfcrenees  94  and  9!>  are  stated  in  terms  of  the  stopping  cross  section 


S (v  ) = - 
1 n 


1 dIl 


i a(Kv1 


dx 


"2V1 


At 


(C-l) 


where  the  neutral  density  n?  is  in  atoms/cm  . The  observed  stopping  cross 
section  includes  a contribution  due  to  elastic  scattering  ("nuclear  stopping"); 
this  is  calculated  to  he  less  than  10  percent  of  the  total  for  speeds  Vj 
greater  than  5 x 10  cm/ sec  (Reference  94).  The  remainder  is  due  to  in- 
elastic scattering  ("electronic  stopping")  and  can  lie  approximated  by  the 
empirical  relation 


S. 


inelv  l 


(v.)  = k v 


o 1 


(C-2) 


where  the  parameters  k^  and  l*  arc  given  in  Table  C-l. 


C-l 


Table  C-1.  Parameters  of  the  stopping  cross  section  due  to  inelastic 
scattering  as  given  by  Sine](vi)  = sinel  is  in 

erg-cm2/atom  and  v,  is  in  cm/sec. 


Colliding 

pair 

Coefficient  k0 

Exponent 

P 

Speed  range, 
cm/sec 

Reference 

N on  air 

4.9*10" 34 

1 .0 

5*107<v1<2.6*108 

94  to  97 

0 on  air 

5. 0*10“ 34 

1.0 

5*107<v1<2.5*108 

94.  95 

Ac  on  air 

3.8*10'35 

1.14 

1 .2*108<v1<3.7*I08 

95,97 

lie  on  lie 

6.1  * 10"36 

1.00 

2.7x108<v1 <4.4*108 

96,  97 

U on  N 

1.8  * 10"37 

l 47 

1 .5*108<v1<1 .0*109 

98 

How  is  S.  , (v.l  related  to  the  momentum- transfer  cross  section 
i no  1 1 

q(v,e)?  We  first  note  that  inelastic  collisions  involve  many  different 
kinetic  energy  losses  c,  depending  on  the  final  states  of  the  colliding 
particles.  We  denote  the  pair  of  final  states  by  i and  treat  the  energy 
losses  as  discrete:  c^,  l ! 1,2,...  . The  rates  of  energy  and  momentum 

transfer  depend  on  two  sums  of  cross  sections: 

l‘inel(v')  - L'  (C"3) 


and  the  quantity  y^e.^  (v.t^) . The  electronic  stopping  cross  section 

S.  . (v,)  can  be  related  to  these  by  Uq.  2-37  with  T = 0,  T = 0,  and 
^mel  1'  1 12 

= 0 and  F.qs.  2- 33c  and  2- 35c: 

wl~ra2  > m? 

WV  * : r?  * ovSuT  W.'V't'- 
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In  evaluate  q.^^jlv)  wo  make  two  assumptions.  At  these  high 
speeds  the  kinetic  energy  tends  to  dissipate  in  many  collisions  with  small 
losses  rather  than  in  a lew  with  hij;  losses,  so  we  can  assume  that 

mm, 

I'd  ~ 2 hynj)  v"  >C-5) 

whenever  t.o  iv,t{)  is  important.  We  further  assume  that  scattering  in 
these  inelastic  collisions  is  in  the  forward  direction, 


t)  <<  l. 


from  the  definitions  of  Mjnej(v)  anil  o^tv.u^)  (f.qs,  C-3,  2-23,  and  2-22) 
one  then  finds  that 

'Y’o 

2JcK.°t (v »c*.)  = (mTImT)-  v qinel(V)  * ,,f  ; tc:'7-) 

t 1 l 

Approximat  ions  C-.r>  and  C •<>  imply  that  the  correction  term  (.indicated  by 

dots)  is  small  because  it  involves  higher  orders  of  the  small  quantities 

■> 

2e^/pV  and  0.  !t  can  be  shown  independently  of  approximations  C-5  and 
C-o  that  this  correction  term  is  negative,  as  indicated. 

Combining  this  result  with  f.qs.  C-4  and  one  finds  that 

lm1+m2)  S.noltv) 


‘‘inellv) 


(ii^rm-,)  k v*’ 


£WV’C*>  = Sineltv)  - 


where  the  neglected  terms  have  the  signs  indicated. 

figures  4-2,  4-3  and  4-5  show  Mjneflv)  for  N,  0,  At,  and 
U ions  on  N calculated  using  liq.  C-S  and  Table  ('-l.  for  ion  molecule 
collisions  we  have  added  the  appropriate  ion-atom  cross  sections  (figures 
i-o  to  4-8). 
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